US  ARMY 
MATERIEL 
rr*MMAND 


TECHNICAL  REPORT  BRL-TR-2736 


•Original  contains  color  v. 
platasx  All  SIic  reproduot- 
lorts  *iu  b»  la  black  and  - 
vhita* 


A  COLOR/GRADIENT  PRE-  AND  POST¬ 
PROCESSOR  FOR  2-D  FINITE  ELEMENT 
STRESS  ANALYSIS  PROGRAMS 


James  II.  Bender 


June  1986 


i 


DTIC 

SLECTEI 
SEP  2  3  9869 


APPROVED  FOR  PURUC  RELEASE;  DtSTRMUTKM  UNLIMITED. 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


l>.>? 

I  V  ’ 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dele  Entered) 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO.j  5  RECIPIENT'S  CATALOG  NUMBER 


Technical  Report  BRL-TR-2736 


4.  TITLE  (end  Subtitle) 

A  COLOR/GRADIENT  PRE-  AND  POST-PROCESSOR  FOR  2-D  j 
FINITE  ELEMENT  STRESS  ANALYSIS  PROGRAMS 


TYPE  OF  REPORT  &  PERIOD  COVEREO 


FINAL  REPORT _ 

S  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORf.; 


B.  CONTRACT  OR  GRANT  NUMBER'*) 


JAMES  M.  BENDER 


9.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-IB 

Aberdeen  Fr ovine-  Ground,  IT)  2100^-5066 


t0.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  »  WORK  UNIT  NUMBERS 

1L162 618AH80 
1L161102AH43 


II.  CONTROLLING  OFFiCP  NAME  ANO  ADDRESS  12.  REPORT  DATE 

U.S.  Army  Ballistic  Research  Laboratory  _ June  1986 _ _ 

ATTN:  SLCBR-DD-T  13  number  of  pages 

Aberdeen  Proving  Ground.  MD  21005-5066 _ 68 _ 

4.  MONITORING  AGENCY  NAME  6  AODRESS</f  different  from  Controlling  Office)  IS  SECURtTYCLASs7fo7fh777»^«j 

UNCLASSIFIED 

IS*,  decl  assific  atPon  downgrading 

SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  I hie  Report) 


Approved  for  Public  Release;  Distribution  Unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  ebetrecl  entered  In  Bloc*  20,  If  different  from  Report) 


11  KEY  WORDS  (Contlnua  on  ravaraa  a/da  II  nacaaaary  and  Idantlfy  by  block  numbar) 

Finite  Element  Color  Code 

Artillery  Projectile  Stress  Analysis 

Pre-processor 
Post-processor 


20.  ABSTRACT  (Continue  «n  reeerem  ride  ft  noceeeery  end  Identify  by  block  number)  j m>C 

A  pre-  and  post-processor  of  finite  element  data  for  use  on  a  color 
graphics  microcomputer  has  recently  been  developed .within  the  Mechanics  & 
Structures  Branch  at  the  U.S.  Army*s  Ballistic  Research  Laboratory.  >The  entire 
data  file,  which  includes  information  for  gridding,  material  properties  and 
loading  schemes,  can  be  generated  and  debugged  locally  on  a  microcomputer  before 
being  sent  to  the  mainframe  computer.  Before  the  emergence  of  advanced,  large 

(continued 


f  tdmnitfy  by  block  numbar) 


ull 


FORM 
I  JAN  73 


EDITION  OF  I  MOV  63  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  TMlS  PA  .E  (When  Dele  Fnlrrnrfl 


SECURITY  CLASSIFICATION  OF  This  PAGEfWIien  Data  Bn  land) 


(continued) 


20.  memory  microcomputers,  the  preprocessing  was  performed  interactively  with 
the  mainframe  computer  which  consumed  a  considerable  amount  of  wall-clock  time. 
Another  advantage  of  pre-processing  on  microcomputers  is  that  the  mainframe  is 
less  burdened  by  interactive  usage  and  is  returned  to  the  batch  mode  allowing 
it  to  perform  more  jobs  in  less  time. 

The  post-processor  transforms  the  output  data  into  a  representative  color, 
denoting  a  component,  and  intensity  of  that  color,  denoting  a  normalized  stress 
level  in  the  element.  This  allows  the  analyst  to  examine  the  state  of  stress 
in  a  multi-component  structure  in  a  single  picture.  For  effective  stress  the 
stresses  are  normalized  with  respect  to  the  material  yield  strength  in  each 
component.  For  the  three  unidirectional  and  shear  stresses  the  elements 
comprising  each  component  are  searched  for  the  maximum  value  and  stresses  are 
normalized  to  that  value.  The  maximum  stress  levels  are  listed  in  a  key  at  the 
top  of  the  display.  This  process  reduces  the  need  for  printed  output  to  a 
backup-only  basis . 


r. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  This  P  AGEfWTien  Data  Entered) 


TABLE  OF  CONTENTS 


Page 


LIST  OF  ILLUSTRATIONS .  5 

I.  INTRODUCTION .  7 

II.  BACKGROUND .  7 

III.  PRE-PROCESSOR .  7 

IV.  POST-PROCESSOR . 8 

V.  CONCLUSIONS . 15 


REFERENCES . 

APPENDIX  A  PREPPY  USER'S  GUIDE 
APPENDIX  B  PROGRAM  LISTING.... 
DISTRIBUTION  LIST . 


19 

21 

27 

63 


LIST  OF  ILLUSTRATIONS 


Figure 


Color-Coded  Finite  Element  Grid  of  the  Rear  Region  of  an 
Artillery  Projectile. .  9 


Numbered  Element  Grid . . . . .  10 


Drawing  of  Projectile  and  Finite  Element  Grid  Reflection....  11 


Loads  on  the  Projectile  During  Launch . .  12 


Compressive  Axial  Stresses, 


Tensile  Radial  Stresses, 


Tensile  Hoop  Stresses . . .  17 


Positive  Shear  Stresses . . . . .  17 


Effective  (Von  Mises)  Stress. 


J 


Accesion  For 


"TIS  CRA&J 
DTiC  TAB 

Unannounced 

Justification 


By . 

Distribution  / 


Availability  Codes 


Dist  !  Avaj'  and/or 
I  special 


'  h*  IN,  m".*  fcN,  •  *  -“**  L.N 


I.  INTRODUCTION 


The  use  of  the  finite  element  method  of  stress  analysis  to  study 
artillery  projectiles  has  been  prevalent  for  at  least  the  past  decade.  The 
complexity  of  these  projectiles  is  growing  steadily  and  likewise  is  the  need 
for  highly  accurate  and  realistic  finite  element  methods.  What  is  presented 
here  is  not  a  new  finite  element  code  but  a  new,  highly  understandable  method 
of  grid  generation  (the  pre-processor)  and  stress  output  display  (the  post¬ 
processor).  The  pre-  and  post-processor  PREPPY  was  written  within  the  U.S. 

Army's  Ballistic  Research  Laboratory  for  a  2-D  axisymmetric  finite  element 

code.1  The  code  takes  full  advantage  of  the  color  graphics  capabilities  of  a 

microcomputer  system. 


II.  BACKGROUND 

Until  the  onset  of  high  technology  warfare  most  artillery  projectiles 
consisted  of  a  high  explosive  housed  in  a  steel  shell  and  capped  with  a 
fuze.  Today's  modern  arsenal  still  includes  high  explosive  types  but  also 
includes  many  cargo-carrying  varieties  in  which  the  payload  is  disseminated 
through  an  expellable  base  at  the  rear  of  the  projectile.  An  ejection  charge 
housed  in  the  nose  cone  is  initiated  by  one  of  many  types  of  fuzes.  The 
expanding  gases  push  on  the  payload  which  in  turn  pushes  on  the  base  and 
shears  the  screw  threads  coupling  it  to  the  body.  The  cargo,  which  could  be 
grenades,  land  mines  or  smoke  pellets,  for  example,  then  fall  to  the  ground, 
so  it  is  obvious  the  many  components  of  these  projectiles  must  be  accurately 
modeled  in  the  finite  element  grid. 

III.  PRE-PROCESSOR 

One  can  imagine  that  a  structure  made  of  three  different  materials,  say, 
steel,  plastic  and  aluminum  can  have  up  to  ten  materials  in  the  finite  element 
domain.  These  include  the  three  real  materials  and  as  many  as  seven 
artificial  materials  which  could  be  tailored  to  behave  as  screw  threads, 
interfaces  or  voids  between  any  two  of  the  materials.  A  natural  course  of 
action  would  be  color  coding  as  in  the  pre-processor  of  PREPPY. 

PREPPY  is  a  computer  program  which  can  be  used  on  a  microcomputer  to  pre- 
process  the  input  data  file  and  send  it  on  its  way  to  the  mainframe  computer 
undisturbed.  The  program  enables  the  user  to  generate  and  de-bug  the  grid 
locally  which  is  less  costly  and  less  time  consuming  than  the  old  method  of 
interacting  with  the  mainframe  computer.  The  mainframe  computer  also  benefits 
in  that  it  is  not  burdened  with  interactive  usage  and  thereby  is  free  to  work 
in  the  batch  mode  for  which  it  was  designed. 


.Tones,  R.M. .  Crouse.  J.G.,  "SAAS  II  Finite  Element  Stress  Analysis  of 
Axisymmetric  Solids  with  Orthotropic,  Temperature  Dependent  Material 
Properties,"  Aerospace  Corporation,  San  Bernadino,  CA,  September  1968. 

Bender,  .T.M. ,  Bums,  B.P..  "Use  of  the  Finite  Element  Method  of  Stress 
Analysis  as  an  Aid  in  the  Design  of  Large  and  Small  Caliber  Projectiles," 
Proceedings  of  the  2nd  Annual  ASME  Computer  Engineering  Conference,  1982. 
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PREPPY  reads  the  input  data  file  in  the  format  identical  to  that  used  in 
the  main  program.  It  creates  a  1000  x  11  array  containing  11  bits  of 
information  about  each  element,  including  the  element  number,  material 
identification  number  (1  through  10),  the  x-  and  y-  coordinates  of  the  four 
corners,  and  the  yield  strength  of  the  element  material.  This  array  is  then 
permanently  stored  on  disk  and  is  available  to  create  color  coded  grids  as  in 
Figure  1  or  numbered  element  grids  as  in  Figure  2.  The  color  coded  grid  is 
also  accessed  by  the  post-processor  for  displaying  output  data. 

Figure  3  illustrates  how  the  manufacturer's  drawings  are  transformed  into 
a  finite  element  grid.  Typically  only  thu  rear  portion  of  projectiles  are 
analyzed  since  that  is  where  body  forces  accumulate  in  the  severe  in-bore 
phase  of  the  launch.  Generally,  the  loads  on  projectiles  consist  of 
propelling  charge  pressure,  axial  acceleration  and  rotating  band  pressure  as 
shown  in  Figure  4.  In  the  case  shown,  the  rotating  bands  exert  an  inward¬ 
acting  radial  pressure  of  about  30  kpsi.  Propellant  gasses  exert  a  pressure 
of  about  50  kpsi  to  the  base  up  to  the  rotating  band.  The  resulting  axial  in¬ 
bore  acceleration  is  approximately  15,000  g's  . 

When  the  user  is  satisifed  that  the  grid  and  loading  scheme  are  in  order 
the  microcomputer  can  be  placed  in  the  terminal  mode  and  transfer  the  file  to 
the  mainframe  computer  for  execution. 

IV.  POST-PROCESSOR 

The  main  program  to  which  the  data  file  is  sent  includes  a  routine  to 
write  selected  output  data  to  a  file  which  is  easily  accessed  by  the 
microcomputer  acting  in  the  terminal  mode.  After  the  file  is  transferred  to 
disk  the  microcomputer  is  switched  back  to  command  mode  where  the  post¬ 
processor  of  PREPPY  will  go  to  work. 

Up  to  this  point  post-processing  was  performed  by  hand  numbering  the 
finite  element  grid  and  searching  through  reams  of  paper,  matching  printed 
data  to  the  grid.  If  one  was  so  inclined,  colored  highlighters  could  be  used 
to  color  in  elements  according  to  material  and/or  stress  level.  This 
procedure  can  be  very  time  consuming  and  tedious,  especially  when  performing 
parametric  studies  which  could  result  in  as  many  as  twenty  output  files.  The 
post-processor  of  PREPPY  has  made  this  method  obsolete.  PREPPY  will 
automatically  generate  colored  stress  output  by  the  following  method.  First, 
PREPPY  accesses  the  element  information  array  developed  during  pre¬ 
processing.  The  eleventh  column  of  this  array  contains  the  color  code  number 
(0  through  1)  for  each  element.  As  examples,  0.1  will  decode  to  the  color 
red,  0.2  decodes  to  orange,  0.5  decodes  to  blue  and  so  forth.  However,  this 
number  (hue)  is  only  one  of  three  numbers  needed  for  full  color  description. 
The  two  remaining  numbers  (also  0  through  1)  are  the  intensity  and 
luminosity.  The  luminosity,  or  brightness,  is  left  at  the  value  1.0  for 
maximum  brightness.  However,  the  intensity  value  will  be  varied  according  to 
the  stress  level.  The  element  will  retain  its  color  code  but  range  in 
appearance  from  white  through  pink  to  red  to  denote  low  through  medium  to  high 
stress  levels  respectively  in  elements  color  coded  by  red.  Red  for  example, 
denotes  aluminum.  In  a  like  manner,  the  other  nine  distinct  colors 
(materials)  are  coded.  The  level  of  intensity  is  determined  by  accessing  the 
1000  x  5  stress  output  array  which  is  transferred  to  the  microcomputer's  disk 
by  the  mainfrime  computer.  This  array  contains  the  axial,  radial,  hoop,  shear 
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Figure  3.  Drawing  of  Projectile  and  Finite  Eleaent  Grid  Reflection 


and  effective  stressJ>  levels  for  each  element.  For  the  effective  stress  the 
value  is  picked  off  the  output  array  and  divided  by  the  yield  point  of  the 
material  of  the  same  element  on  the  element  information  array.  This  quotient, 
say  0.7  for  example,  will  create  an  element  with  70%  intensity  of  the 
designated  color.  Tables  1  and  2  are  partial  lists  of  the  element  information 
array  and  the  stress  output  array  respectively,  and  illustrate  the  procedure. 


TABLE  1 


Element 

Radial 

Hoop 

Axial 

Shear 

Ef  fective 

# 

Stress 

Stress 

Stress 

Stress 

Stress 

( psi) 

(psi) 

( psi) 

( psi) 

(psi) 

345 

36,264 

16,035 

107  ,892 

16,514 

121  ,031 

346 

38,392 

14,901 

106,041 

16,617 

107  ,389 

347 

40,392 

138,771 

101 ,581 

17,722 

64,993 

TABLE 

2 

Element 

I.D. 

Coord i nates 

Yield 

it 

it 

LL 

UL 

UR 

LR 

Stress 

(psi) 

345 

2 

x»y 

*,y 

x.y 

x,y 

140,000 

346 

2 

*,y 

x,Y 

x,y 

x,y 

140,000 

347 

4 

*,y 

*,y 

x.y 

x,y 

65,000 

Referring  to 

Tables 

1  and 

2,  for  element 

346  the 

intensity 

would  be 

107,389/140, 

000  or  0 

.767  ■ 

which  decodes  to 

orange 

at  about 

77%  color 

intensity.  Element  347  decodes  to  blue  at  64,993/65,000  or  nearly  100%  color 
intensity.  This  element's  stress  is  near  the  yield  point.  Should  the  value 
exceed  102%,  an  "x"  is  placed  on  the  element,  indicating  that  perhaps  an 
elastic-plastic  computation  is  warranted. 

For  stress  in  a  uniaxial  direction  (axial,  radial,  hoop,  or  shear)  the 
elements  must  be  grouped  according  to  ident i f i cat i on  number.  Each  group  is 
searched  for  the  maximum  value  in  each  stress  direction.  These  maximum  values 
are  listed  in  a  key  at  the  top  of  the  stress  plot.  The  same  procedure  is  used 
as  in  effective  stress  to  compute  intensity  except  that  the  group  maximum 


Hill,  R. ,  The  Mathematics1  Theory  of  Plasticity,  Oxford  University  1'ress, 
London,  England,  1950. 

Thomas,  T.Y. ,  Plastic  Flow  and  Fracture  in  Solids, 

York,  NY,  1961. 
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value  replaces  the  yield  stress  in  the  computation.  Tables  3  and  4 
illustrate  the  procedure  where  Table  3  contains  the  group  of  elements  with 
identification  number  6  and  Table  4  contains  pertinent  portions  of  the  element 
information  array. 

TABLE  3. 


Element 

Radial 

Hoop 

Axial 

Shear 

Effective 

# 

Stress 

Stress 

Stress 

Stress 

Stress 

(psi) 

(psi) 

(psi) 

(psi) 

(psi) 

23 

21,962 

16,123 

(60.0071 

1692 

63,222 

24 

23,091 

16,242 

59,066 

2023 

63,359 

25 

27,211 

16,397 

57,243 

2394 

63,494 

196 

29,997 

17,291 

54,129 

2711 

64,093 

197 

31 ,035 

18,121 

52,449 

2638 

64,125 

198 

(33,323| 

(l9.892| 

50,030 

2469 

64,210 

TABLE  4 


Element 

i.n. 

Coordinates 

Yield 

if 

if 

LL 

UL 

UR 

LR 

Stress 

(psi) 

23 

6 

x,y 

x.y 

x.y 

x.y 

65,000 

24 

6 

x.y 

x.y 

x.y 

x.y 

65,000 

25 

6 

x.y 

x.y 

x.y 

x.y 

65,000 

196 

6 

x.y 

x.y 

x.y 

x.y 

65,000 

197 

6 

x.y 

x.y 

x.y 

x.y 

65,000 

In  Table  3,  the  maximum  values  are  boxed  in.  Element  23  would  be  green  at 
21,962/33,323  or  66%  intensity  for  the  radial  stress  plot.  If  the  axial 
stress  plot  were  generated  element  23  would  be  green  at  60,007/60,007  or  100% 
intensity  since  this  is  the  highest  stressed  element  in  the  axial  direction. 
Table  4  is  used  to  obtain  the  coordinates  of  the  particular  element. 

Thus  far  only  positive  or  tensile  stresses  have  been  examined.  However, 
the  MATRIX  SEARCH  routine  includes  provisions  to  search  for  minimum  values 
(maximum  compression)  as  well  and  groups  these  in  a  like  manner  as  the  tensile 
stresses.  .Just  as  .negative  (compression)  values  are  ignored  during  the  search 
for  maximum  tension  in  each  group,  positive  values  are  ignored  during  the 
search  for  maximum  compressive  stresses.  On  the  actual  plots,  the  ignored 
elements  are  colored  in  at  full  intensity  but  at  20%  luminosity  for 


completeness  of  the  figure.  Stressed  elements  appear  much  brighter.  Figures 
5  through  9  are  examples  of  PREPPY's  stress  plots. 


V.  CONCLUSIONS 

It  can  be  seen  that  there  are  nine  possible  stress  plots  for  each 
computer  run:  radial,  hoop,  axial,  and  shear  all  in  either  compression  or 
tension  (shear  merely  shows  direction)  and  the  effective  stress.  In  most 
cases  tensile  stresses  are  given  more  attention  since  they  tend  to  impact 
crack  propagation  and  critical  defect  size.^  In  practice  over  the  past  year 
PREPPY  has  been  found  to  help  the  seasoned  engineer  by  making  efficient  use  of 
time  and  as  a  visual  aid  in  communications  with  the  higher  echelons.  PREPPY 
has  also  proven  itself  as  a  valuable  training  tool  for  new  engineers  to  help 
them  visualize  and  recognize  stress  patterns  in  structures. 
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Hertzberg,  R.W. ,  Deformation  and  Fracture  Mechanics  of  Engineering 
Materials,  John  Wiley  &  Sons,  NY,  1976. 
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APPENDIX  A 
PREPPY  USER'S  GUIDE 


USER'S  GUIDE  FOR  PREPPY 


1.  CREATE  A  DATA  FILE  JUST  AS  IF  YOU  WERE  RUNNING  ON  THE  CDC.  THIS  CAN 
BE  DONE  USING  RUFUS.  THE  TWO  TYPES  OF  RUN  STREAMS  ARE  LISTED  IN  THE 
APPENDIX.  EACH  TYPE  HAS  A  SPEC  I F f  C  FUNCTION  WHICH  WILL  BE  DISCUSSED 
LATER.  DO  NOT  INCLUDE  THE  END-OF-RECORD  MARKS. 

2.  AFTER  THE  DATA  FILE  IS  CREATED  DO  A  LOAD"PREPPY: D" .  THE  FIRST  TIME 
THROUGH  YOUR  MENU  SELECTION  WILL  BE  “D“  FOR  DATA  FILE.  THIS  ROUTINE 
WILL  COMPUTE  AND  GENERATE  YOUR  GRID.  IF  YOU  ARE  HAPPY  WITH  THE  GRID 
YOU  CAN  STORE  IT  AS  A  “GRID"  FILE.  ACTUALLY,  YOU  ARE  PLOTTING  INTO 
AN  ARRAY  WHICH  WHEN  RECALLED  TAKES  A  FRACTION  OF  THE  TIME  TO  REPLOT 
COMPARED  TO  GENERATING  THE  GRID.  YOUR  GRID  WILL  APPEAR  ON  THE  CRT, 
HOWEVER  YOU  ARE  GIVEN  THE  OPTION  OF  A  HARD  COPY.  THE  HARD  COPY  WILL 
BE  A  LARGE  PLOT  WHICH  FITS  ON  USED  COMPUTER  PAPER  OR  LARGE  HP  PAPER 
RATHER  WELL.  PLEASE  BEGIN  THE  NAME  OF  THE  GRID  FILE  WITH  A  "G"  SUCH 
AS  “RING"  BECOMES  "GRING"  AS  A  GRID  FILE. 

3.  AFTER  EXERCISING  YOUR  GRID  STORAGE  OPTION,  PREPPY  WILL  CONTINUE  AND 
BEGIN  TO  FILL  IN  THE  ELEMENTS  WITH  A  COLOR  ACCORDING  TO  YOUR  MATER¬ 
IAL  BLOCK  ASSIGNMENT  CARDS.  KEEP  IN  MIND  THAT  GREEN  (2),  BLUE  (4), 
AND  RED  (6)  PHOTOGRAPH  THE  BEST.  AS  AN  ATTEMPT  TO  STANDARDIZE  THE 
COLOR  CODING  SCHEME,  TRY  TO  RESERVE  BLUE,  FOR  STEEL  AND  RED,  FOR 
ALUMINUM.  YELLOW  (1)  DOES  NOT  SHOW  UP  WELL  SO  USE  IT  FOR  AIR.  SAAS 
DOES  NOT  CARE  WHAT  ORDER  THE  MATERIAL  BLOCK  ASSIGNMENT  CARDS  OR  THE 
MATERIAL  INFORMATION  CARDS  ARE  IN.  HOWEVER,  UNLIKE  RUNNING  SAAS  ON 
THE  CDC  PREPPY  IS  PARTICULAR  ON  ASSIGNING  MATERIAL  TO  NON-EXSISTENT 
ELEMENTS,  SO  COLOR  IN  ONLY  THE  ELEMENTS  YOU  HAVE.  GRANTED,  THIS 
CREATES  A  FEW  MORE  LINES,  BUT  THE  OVERALL  TIME  SAVINGS  IS  SIGNIFI¬ 
CANT  AND  MAKES  NO  DIFFERENCE  TO  SAAS  ON  THE  CDC. 

4.  IF  YOU  ARE  HAPPY  WITH  THE  COLOR  CODED  GRID,  PREPPY  WILL  ASK  YOU  IF 
YOU  WANT  TO  STORE  A  SUPERGRID.  THE  SUPER  GRID  FILE  CONTAINS  ALL  THE 
INFORMATION  NECESSARY  FOR  RECREATING  THE  COLOR-CODED  GRID  QUICKLY, 
THE  NUMBERED  ELEMENT  GRID,  AND  THE  STRESS  PLOTS.  IT  IS  SUGGESTED 
THAT  THE  NAME  OF  THE  SUPERGRID  FILE  START  WITH  "S"  SO  THAT  IT  IS 
EASILY  RECOGNIZED  AS  SUCH  (EXAMPLE:  THE  SUPERGRID  FILE  OF  "GAGE" 
BECOMES  "SGAGE").  KEEP  IN  MIND  THAT  THE  9845  WILL  SPEND  ABOUT  FOUR 
MINUTES  COMPUTING  COORDINATES  FOR  300  ELEMENT  GRID. 

5.  TO  OBTAIN  THE  GRID  FILE  SELECT  "G"  FROM  THE  PREPPY  MENU.  ENTER  THE 
NAME  OF  THE  GRID  FILE.  SELECT  THE  SCALE:  "FULL"  MEANS  THE  ENTIRE 
FIGURE,  "1/2"  MEANS  THE  FIRST  HALF,  "2/4"  MEANS  THE  SECOND  QUARTER 
ETC,  ETC.  IF  YOU  WANT  A  HARD  COPY  SELECT  THE  SIZE  OF  THE  PAPER  AND 
THE  SCALE.  IF  NOT,  THE  GRID  WILL  APPEAR  ON  THE  SCREEN. 

6.  IF  YOU  WANT  A  NUMBERED-ELEMENT  GRID,  YOU  MUST  HAVE  FIRST  CREATED  A 
SUPERGRID  FILE.  SIMPLY  ENTER  THE  NAME  OF  THE  SUPERGRID  FILE,  THEN 
SELECT  THE  SCALE.  PLACE  AN  OVERSIZE  SHEET  OF  PAPER  (COMPUTER)  ON 
THE  PLATEN.  THIS  PROGRAM  WILL  NOT  RUN  UNLESS  THE  PLOTTER  IS  TURNED 
ON.  BECAUSE  OF  JAGGIES  THE  NUMERALS  DO  NOT  SHOW  UP  WELL  ON  THE  CRT. 

7.  IF  YOU  ARE  NOW  READY  FOR  STRESS  PLOTS,  YOU  MUST  HAVE  FIRST  USED  THE 
LONG  RUNSTREAM  (APPENDIX).  THE  THIRD  ENTRY  ON  LINE  14  PROVIDES  FOR 
THE  OUTPUT  FILE  USED  BY  PREPPY  TO  GENERATE  THE  STRESS  PLOTS.  NOW  IS 
THE  TIME  TO  TRANSFER  THE  DATA  FILE  TO  THE  CDC.  DON'T  FORGET  TO  ADD 
IN  THE  END-OF-RECORD  MARKS  BEFORE  THE  HEADER  LINE  AND  THE  END.  WHEN 
THE  PROGRAM  IS  EXECUTED  THE  OUTPUT  FILE  (USUALLY  STARTED  WITH  A  P) 
WILL  AUTOMATICALLY  BE  CREATED.  THIS  FILE  MUST  BE  TRANSFERRED  OVER 
TO  THE  9845'S  DISK.  SELECT  "S"  FROM  THE  PREPPY  MENU  AND  FOLLOW  THE 
INSTRUCTIONS.  MAKE  ABSOLUETIY  CERTAIN  THAT  THE  OUTPUT  OR  P-FILE  IS 
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COMPATIBLE  WITH  THE  SUPERGRID  FILE.  TO  HELP  YOU  KEEP  FI 
USE  THE  PREPPY  LOG  BOOK  IN  RM  107,  OR  CREATE  YOUR  OWN 
CAN  CREATE  UP  TO  NINE  STRESS  PLOTS  FOR  EACH  RUN:  RADIAL 
OR  HOOP,  AXIAL  SHEAR  (ALL  IN  TENSIONS)  OR  COMPRESSION< 
VON  MISES  OR  EFFECTIVE  STRESS.  FOR  THE  UNIAXIAL  DIRECTI 
EST  STRESSED  ELEMENT  IS  THE  DARKEST  COLORED  ONE  IN  THE 
ALL  OTHER  STRESSES  ARE  NORMALIZED  TO  THAT  VALUE.  THE 
THAT  STRESS  IS  LISTED  IN  A  KEY  AT  THE  TOP  OF  THE  PLOT. 
ECTIVE  STRESS  ALL  STRESSES  ARE  NORMALIZED  WITH  RESPECT 
POINT  OF  THE  PARTICULAR  MATERIAL  INVOLVED.  THE  DARKEST 
MENTS  ARE  THOSE  WHOSE  VON  MISES  STRESS  IS  VERY  NEAR  THE 
OF  THE  MATERIAL.  SHOULD  AN  “X"  APPEAR  ON  THE  ELEMENT  I 
STRESS  HAS  EXCEEDED  THE  YIELD  STRESS  BY  TWO  PER  CENT. 


LES  STRAIGHT 
SYSTEM.  YOU 
,  TANGENTIAL 
-),  AND  THE 
ONS  THE  HIGH 
COLOR  GROUP. 
MAGNITUDE  OF 
FOR  THE  EFF- 
TO  THE  YIELD 
COLORED  ELE- 
Y I  ELD  POINT 
TS  EFFECTIVE 


8.  SHOULD  YOU  NEED  A  BLACK-AND-WHITE  VERSION  OF  (7),  SELECT  "BW"  FROM 
FROM  THE  MENU.  ONCE  THE  PLOT  IS  COMPLETED  YOU  CAN  DUMP  THE  SCREEN 
ON  TO  THE  DOT  MATRIX  PRINTER  BY  EXECUTING  A  "DUMP  GRAPHICS".  DON'T 
FORGET  TO  LABEL  THE  COMPONENTS  OF  THE  GRID  SINCE  YOU  NO  LONGER  HAVE 
THE  COLOR-CODING  ADVANTAGE. 


PROCEDURE  FOR  CDC< - >9845  TRANSFERS 


1.  LOAD  "CYBER:D",1  (OR  HIT:  k 3  FROM  ANPACK,  THEN  SELECT:  2,  THEN 
ENTER:  PAX  THEN  MFA  OR  JUST  MFA > 

2.  TURN  ON  BLUE  BOX. 

3.  HIT:  CONT. 

4.  LOG  IN. 

5.  HIT:  k 1 3 . 

6.  HIT:  STORE 

7.  CHANGE  "TEST:T14"  TO  f i 1 ename : C 1 2 

8.  CHANGE  "SIZE=010"  TO  “SIZE«300" 

9.  HIT: STORE 

10.  ENTER:  TRMDEF,PW=81 

-  9845  TO  CDC  - 

11  ENTER:  NEW/LFN1  NOTE:  THE  ARROW  KEYS  ONLY  APPEAR 

12.  ENTER:  TEXT  TO  TO  WORK  WHEN  TALKING 

13.  HIT:  SHIFT  k  1 4  WITH  THE  CDC.  YOU  MUST  USE 

14.  HIT:  N  THE  SPACEBAR  OR  BACKSPACE. 

15.  AFTER  UPLOAD  HIT:  CONTROL  T 

16.  HIT:  CONT 

17.  ENTER:  SAVE,LFNl=f ilename  POPULAR  XEDIT  COMMf 


11.  HIT: 

12.  ENTER 

13.  ENTER 

14.  TYPE: 

15.  HIT: 

16.  HIT: 

17.  AFTER 
HIT: 


CDC  TO  9845  - 

SHIFT  k 1 2 
?:  GET  f  i  1  ename 
*:  ASSIGN, TT,B 
COPY, f i 1 ename, B 
SHIFT  k 1 5 
CONT 

*  FILE  IS  READ  IN 
SHIFT  k 1 5 


t: 

b: 

L^str i ngs: 


DR: 

C 'st  r i ngl ' st r i ng2/n: 
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POPULAR  XEDIT  COMMANDS 

PRINTS  n  LINES 
MOVES  POINTER  n  LINES 
DELETES  n  LINES  STARTING  AT 
CURRENT  POINTER  POSITION 
MOVES  POINTER  TO  TOP  OF  FILE 
MOVES  POINTER  TO  BOTTOM  OF  FILE 
LOCATES  FIRST  OCCURANCE  OF 
string 

WRITES  EOR  MARK  ABOVE  CURRENT 
POINTER  POSITION 
DELETES  THE  EOR  MARK 
CHANGES  stringl  TO  string2 
n  TIMES  (1  IS  DEFAULT) 


1  4_*  AJI  a*  ^  a  \  •  a  1. 


APPENDIX 


SHORT  RUNSTREAM: 

JAMES, STMFZ, T800,P5. 

ACCOUNT, account . 

ATTACH, SAASII, ID=KOKINAKI 

BEGIN, RUNS A AS, SAASII,PL=150000,SHORT=YES. 

JIM  BENDER  B390  ROOM  107  X6116 


SAMPLE 

FINITE 

ELEMENT 

GRID 

100 

-2  9 

90 

0  72 

15768 

5 

9 

8. 

.05 

SAMPLE 

FINITE 

ELEMENT 

GRID 

22 

54  71 

1 

18 

1  1 

0.000 

0.690 

6  1 

1.490 

0.690 

1 

6  1 

1.490 

0.690 

9  1 

2.200 

0.000 

1 

9  1 

2.200 

0.000 

22  1 

2.812 

0  ■  000 

1 

22  1 

2.812 

0.000 

22  8 

2.992 

1 . 253 

1 

22  8 

2.992 

1.253 

22  12 

3.  033 

1 . 538 

1 

LONG  RUNSTREAM*. 


JAMES, STMF2,T900,P5. 

ACCOUNT , account . 

NDFILEC1) 

COPYSP, INPUT, OUTPUT. 

REWIND, INPUT. 

FILE, TAPE 10,RT*2,BT=C,FL=80. 

BEGIN, ATTACH, PLOTLIB. 

ATTACH, LGO, SAAS 1 1 SHORTLGO , I D=KATHY , PW=BRUCE , MR= 1 . 

LGO, PL-100000. 

EXIT, U. 

REWIND, TAPE9. 

REWIND, TAPE10. 

COPY, TAPE9, OUTPUT. 

SAVEPF,TAPE10,Pf i 1 «■ , CH=PD902B , PN  =  PD , TY=REPLACE , UN=name , UP=passwd , ST=MFA . 
BEGIN, PLOT, CALCOMP, TAPE3. 


EXIT. 

BEGIN, PLOT, CALCOMP, TAPE3. 

JIM  BENDER  B390  RM107  X6116 


SAMPLE 

FINITE 

ELEMENT 

GRID 

* 

100 

-2  7 

37 

0 

72. 

32. 

5 

9 

20. 

0.  10 

SAMPLE 

FINITE 

ELEMENT 

GRID 

• 

24 

55  48 

2 

14 

1  1 

.000 

0.  500 

24 

1 

1 . 062 

1 . 350 

24  1 

1.062 

1.350 

24 

1 1 

1.062 

1 . 900 

24  11 

1.062 

1.900 

24 

15 

1.062 

2.650 

24  15 

1.062 

2.650 

24 

30 

1.G62 

5.860 

24  34 

1.062 

6.610 

24 

48 

1 .  @62 

10.350 

32. 


000 


3.426 


4 

1 

1 

1 

1 
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PRINTER  IS  16 
OPTION  BASE  1 


PROGRAM:  PREPPY,  WRITTEN  BY  JAMES  M.  BENDER  MSB,  IBD,  BRL  ON 
NOVEMBER  6,  1981.  THIS  SECTION  OF  THE  PROGRAM  GENERATES  THE 
GRID  FROM  THE  RAW  DATA  FILE  AND  CREATES  THE  GRID  FILE  AND  THE 
SUPERGRID  FILE  FOR  LATER  USE.  USER'S  MANUALS  ARE  AVAILABLE  FOR 
USE  OF  THE  PROGRAM  ON  THE  hp9845  COLOR  MICRO  COMPUTERS. 


REVISION  DATES 


110  ! 

JUL 

19, 

1982 

MAR 

19, 

1984 

120  ! 

SEP 

7, 

1982 

MAY 

6, 

1984 

130  ! 

FEB 

1. 

1983 

MAY 

30, 

1984 

140  ! 

MAY 

16, 

1983 

JUL 

31, 

1984 

150  i 

AUG 

13, 

1933 

MAR 

25, 

1985 

160  ! 

OCT 

25, 

1983 

AUG 

29, 

1985 

170  ! 

JAN 

30, 

1984 

INTEGER  Nc  ode ,  Mat  er  i  al  ,  I  m  i  n ,  I  max , Jm i n , Jmax , Code 

COM  Ar<26,  101>,Az<26,  101>,Ncode<26,  1 0 1 > , E 1 em < 900 ,  ll),Material  <26,  101 ) 
COM  Nmel  (100),Nemax<100),Nemin<100),Npmin<100),Npmax(100) 

COM  Numt c , E J emax, Y i e 1 d< 10) 

COM  Imin<100),  Imax<100),Jmin(26),Jmax<26),Maxi  ,  Max  j  ,  Nmt  1  ,Nbc 
COM  Numnp , Nume 1 , Nummat 

COM  R< 1000) , Code < 1000 > , Xr < 1 000 ) , Z < 1 000 ) , Xz< 1000) 

COM  Kp,Rp<250),2p<250) 

COM  Ne 1 e ( 26 ) 

COM  1 1 ou, Ihi gh, J1 ou, Jhi gh 
COM  Npp,Nseg 

COM  Fil e$C 15] , Li ne*C81  )  , Que*, P* 

DIM  L i ne$  < 200 ) C 8 1  ] 

PRINT  PAGE 

PRINT  "  Welcome  to  PREPPY.  Please  select,  from  the  menu  below." 


PRINT  "  Welcome  to  PREPPY. 
PRINT 

PRINT  "Choice  Meaning 

PRINT 

PRINT  "  s D '  Data  F 

PRINT  " 

PRINT 

PRINT  "  'G'  Grid  F 

PRINT  " 

PRINT 

PRINT  "  'S'  Supergri 

PRINT  " 

PRINT  " 

PRINT 

PRINT  "  SBW  Black  an 

PRINT  "  Plo 

PRINT  " 

PRINT 

PRINT  "  EN'  Element 

PRINT  " 

INPUT  "Enter  Se 1 ec t i on " , Type# 

PRINT  PAGE 


Dat  a  File 


Grid  File 


Supergr id  File 


Black  and  White- 
Plots 


Element  Numbers 


Desc  r i pt i on " 

File  from  which  basic  grid" 
i s  der i oed. " 

Stored  picture  of  the  wire" 
frame  grid. " 

Cont ai ns  t he  i nf or mat i on " 
necessary  for  plotting  color-" 
coded  grids  and  stress  plots." 

Enables  stress  plots  to  be" 
plotted  in  black  and  white  on" 
the  dot  matrix  printer" 

Sequent i a 1 1 y  n umbers  el e m ents " 
on  large  paper  on  the  plotter" 
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550  IF  Typet="G"  THEN  LORD  "GRIL:D“,1 

560  IF  T ypet= " S "  THEN  LORD  "PREPl:D",l 

570  IF  Typet="BW"  THEN  LORD  "PREP2:D",1 

580  IF  Typet="EN"  THEN  LORD  "ELNUM: D"  , 1 

590  INPUT  "Hhich  Data  file  would  you  like  to  access?" , Fil et 
600  PRINT  "You  haue  accessed  ";Fi1et 
610  Fi  1 et=Fi  1  et8c" :  C12" 

620  ASSIGN  #1  TO  Fi 1 et 

630  !  ************************************************************************ 

640  !  *  MESH  CONTROL  INFORMATION  * 

650  !  ************************************************************************ 

660  READ  # 1 ; Li  net 

670  READ  #1 ; L i net 

680  READ  #1 ; Li  re* 

690  Run=16 

700  IF  Li net! 1 , 3]=" ATT"  THEN  Run=3 
710  FOR  1=1  TO  Run 

720  READ  ttl ; Li net  (  I ) 

730  NEXT  I 

740  READ  #1; Li  net 

750  Numtc=VAL(Li netC6, 10]  ) 

760  IF  Numt c  =-2  THEN  Numtc=l 
770  Nummat =VAL <L i net C 1 1 , 1 5] ) 

780  Numpc=VAL<Li nett  16, 20]  ) 

790  Numsc=VAL<Li netC21 , 25]  ) 

800  READ  #1 ; Li  net 

810  READ  #1 ; Li  net 

820  READ  #1 ; Li  net 

830  Maxi =  VAL<Li net  I 1,5]) 

840  Maxj=VAL<LinetC6, 10]) 

850  Nseg=VAL<LinetCll, 15]) 

860  Nbc=VAL<LinetI16,203) 

870  Nmtl=VAL(LinetC21,25]) 

880  IF  Pt»"Z"  THEN  GOTO  930 

890  INPUT  "Do  you  want  an  input  1 i st i ng? < Y/N ) " , Quet 


900 

DISP  " 

*  Reading 

in  Mesh  Data  *" 

910 

IF  Quet= “ Y " 

THEN  PRINTER 

IS  0 

920 

IF  Quet="Y" 

THEN  PRINT  " 

Max  I 

= " ; Max i ; " 

Max  J  =" ; Max j ; " 

#  Line 

Seg 

went 

cards  *";Nseg 

921 

IF  Quet*" Y" 

THEN  PRINT 

922 

IF  Quet="Y" 

THEN  PRINT  " 

1 1 

J1  R 1 

Z1  12  J2  R2 

Z2 

13 

J3 

R3  Z3 

Ipt i on" 

930  !  ************************************************************************ 

940  !  *  INITIALIZE  * 

950  !  ************************************************************************ 

960  Kp= 1 

970  Iseg=-1 

980  Pi =PI 

990  MAT  Ncode=ZER 

1000  MAT  Ar=ZER 

1010  MAT  Az=ZER 

1020  MAT  Jmax=ZER 

1030  MAT  Jmi n=(Max j ) 

1040  MAT  Imin=<Maxi) 

1050  MAT  I max=ZER 
1060  MAT  Mat  er i al =ZER 
1070  MAT  E 1 *m*ZER 

1080  !  ************************************************************************ 
1090  !  LINE  SEGMENT  CARDS 

1100  !  ************************************************************************ 

1110  Iscg=Is*g+l 

1120  IF  I seg=0  THEN  1190 

1130  IF  <  Ipt  i on<  =3 )  OR  <Iption=6)  THEN  GOTO  1190 

1140  Rp<  Kp ) =R3 
1150  Zp<Kp)«Z3 
1160  GOTO  1190 
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1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
.  DDD, 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 


Rp  <  Kp  ) *R2 
Zp<Kp>»Z2 

IF  Iseg=Nseg  THEN  GOTO  Interior 
Ic*3 

I1=I2=I3=J1=J2=J3=R1=R2=R3=Z1=22=Z3=0 
READ  «l;Line* 

Ii=VAL<Line*Cl,33> 

Jl=VAL<Line$[4,63> 

R1 =VAL  <L i ne*[ 7, 143) 

Zl=VAL<Line$[15,223) 

Ipt i on=VAL<Li ne$C67,713) 

IF  lption=0  THEN  Next_ic 
I2=VAL<Line$C23,253) 

J2*VAL<Li ne$C26, 283  > 

R2=VAL<Line*C29,363) 

Z2*VAL(Line*f37,443> 

IF  <Iption<>4>  AND  (Iption<>3>  AND  <Iption<>13>  THEN  Next__ic 
R3=VAL<Line$[51,583) 

Z3=VAL<Li ne*C59, 663  ) 

IF  Iption=4  THEN  Next_ic 
I3»VAL<Line$C45,473> 

J3=VAL<Line*[48,503> 

Next_ic:  Ic«Ic+l 

IF  <Que$="N" )  OR  <P$="Z">  THEN  GOTO  1430 

PRINT  USING  Image;ll,Jl,Rl,Zl,I2,J2,R2,Z2,I3,J3,R3,Z3,Iption 
Image: IMAGE  IX,  DDD,  IX,  DDD, IX, DD. DDD, IX, DD. DDD, IX, DDD, IX, DDD, IX, DD. DDD, IX, DD 
IX, DDD, IX, DDD, IX, DDD. DDD, 1 X , DDD . DDD , 3X , DD 
IF  lption=0  THEN  GOTO  Diagonal 
IF  Ipt i on=l 1  THEN  GOTO  1680 
IF  <Iption=5)  OR  <Iption=l)  THEN  Straight 
IF  Iption=2  THEN  GOTO  Diagonal 
IF  Iption=13  THEN  GOTO  Diagonal 
Rp (  Kp ) =R 1 
Zp<Kp)=Zl 
GOTO  Diagonal 

I 

!  *****  POINTS  ***** 

j 

Rp<Kp>=Rl 
Zp ( Kp ) =Z 1 
Kp=Kp+ 1 

IF  Iseg+l<Nseg  THEN  GOTO  Diagonal 
Rp  <  Kp ) *Rp  < 1 ) 

Zp  <  Kp ) =Zp ( 1  > 

GOTO  Diagonal 

Straight:  !  *****  FOR  IPTION  =  1  ***** 

! 

Rp  <  Kp ) =R 1 
Zp ( Kp ) *Z 1 
Rp ( Kp* 1 ) =R2 
Zp ( Kp* 1 ) =Z2 
Kp*Kp*2 

IF  Iption=ll  THEN  Iption=l 


D i agonal : 


i  *****  FOR  IPTION 


0,2  or  13  ***** 


IF  1 1  *- 1  THEN  GOTO  Interior 

IF  Iption=13  THEN  Ipt=13 

IF  Iption=13  THEN  Iption=3 

Ipt i on=Ipt i on+1 

Ar< 1 1 , J1 )*R1 

Az< 1 1 , J 1 >  *Z 1 

Ncode< 1 1 , J 1 >  =  1 

CALL  Mn i mx  < 1 1 , J 1 ) 

ON  Ipt ion  GOTO  1 1 10, 1840, 1840, Ci rc 1 e, Ci rc 1 e, 1840 


>r 


Msnfi 


1820  !  GENERATE  STRAIGHT  LINES  ON  BOUNDARY 

1830  !  ************************************************************************ 

1840  Di *ABS< 12-11) 

1850  D j=A8S< J2-J 1 > 

1860  Ar<I2,J2)=R2 

1870  Az< 12, J2)=Z2 
1880  Ncode<I2, J2)*l 
1890  CALL  Mni mx< 12, J2) 

1900  Istrt=I 1 

1910  Istp=I2 

1920  JstrfJl 

1930  J*tp=J2 

1940  Diff=MAX<Di ,Dj> 

1950  Iter«Di f f-1 
1960  I i nc=0 
1970  Ji nc=0 

1980  IF  12011  THEN  I  i  nc  =  <  1 2- 1 1  ) /ABS <  1 2- 1 1 ) 

1990  IF  J20J1  THEN  J i  nc  =  ( J2- J 1 ) /ABS <  J2- J 1 ) 

2000  Kappa=l 

2010  IF  <12011)  AND  <J2<>J1)  AND  <Iption<>3)  THEN  Kappa=2 

2020  IF  Kappa=2  THEN  Diff=Diff*2 

2030  Rinc=<R2-Rl)/Diff 

2040  2inc=<22-Zl)/Di ff 

2050  ! 

2060  !  *****  CHECK  FOR  INPUT  ERROR  ***** 

2070  ! 

2080  IF  <Kappa<  >2 )  OR  <Di=Dj>  THEN  GOTO  2130 
2090  GOTO  1110 
2100  ! 

2110  !  *****  INTERPOLATE  ***** 

2120  ! 

2130  1=11 

2140  J=J1 

2150  FOR  M=1  TO  Iter 

2160  IF  <Iter=8>  AND  <Iption=2>  THEN  GOTO  2550 

2170  IF  < lter=0)  AND  <Iption=6>  THEN  GOTO  2550 

2180  IF  Kappa=2  THEN  GOTO  2280 

2190  Iol d=I 

2200  1=1+1 i nc 

2210  Jol d=J 

2220  J=J+J i nc 

2230  Ar<I, J)=Ar<Iold, Jold)+Rinc 
2240  Az< I , J)=Az< Iold,Jold)+Zinc 
2250  CALL  Mni mx< I , J) 

2260  Ncode< I , J)=l 
2270  GOTO  2550 

2280  IF  < 1 1 >12)  AND  <Iption=2>  THEN  GOTO  2430 
2290  IF  ( I 1< 12)  AND  <Iption=6)  THEN  GOTO  2430 
2300  Iol d*I 
2310  I=I+Iinc 

2320  Ar< I , J)=Ar< Iold, J)+R j  nc 
2330  Az<I, J)=Az<Iold, J)+2inc 
2340  Ncode< I , J)=l 
2350  CALL  Mni mx< I , J) 

2360  Jold=J 
2370  J»J+Jinc 

2380  Ar< I , J)*Ar < I , Jol d)+Ri nc 
2390  Az<I, J)=Az<I, Jo1d)+2inc 
2400  Ncode< I , J)=l 
2410  CALL  Mnimx< I , J) 

2420  GOTO  2550 
2438  Jold=J 
2448  J«J+Jinc 

2450  Ar < I , J)«Ar < I , Jol d) +Ri nc 
2460  Az< I , J)=Az< I , Jol d)+Zi nc 
2470  Ncod*< I , J)«l 
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2480  CALL  Mni mx< I , J> 

2490  Io1d*l 
2500  1*1+ line 

2510  Ar<  I , J>*Ar< Io1d,J>+Rinc 
2520  Az< I , J>=Az< I > i d, J)+Zi nc 
2530  NcodeC I , J  >  *  1 
2540  CALL  MnimxC I ,  J) 

2550  NEXT  M 

2560  IF  Kappa*l  THEN  GOTO  1110 

2570  IF  ( 1 1 >I2>  AND  <Iption=2>  THEN  GOTO  2640 

2580  IF  < 1 1 < 1 2 >  AND  <Iption=6>  THEN  GOTO  2640 

2590  Io1d*I 

2600  Isl+Iinc 

2610  Ar<  I , J>=Ar< Io1d,J)+Rinc 
2620  Az<I, J>=Az<Io1d, J>+Zinc 
2630  GOTO  2680 
2640  Jol d=J 
2650  J*J+Jinc 

2660  Ar x I , J)*Ar( I , Jol d)+Ri nc 
2670  AzCI, J>*Az<I, Jold>+Zinc 
2680  NcodeC I , J)-l 
2690  CALL  Mri mx< I , J ) 

2700  GOTO  1110 

2710  !  ************************************************************************ 
2720  !  *  GENERATE  CIRCULAR  ARCS  ON  BOUNDARY  * 

2730  !  ************************************************************************ 
2740  Circle:  ! 

2750  Ar<I2,J2>*R2 

2760  Az<12, J2>*Z2 

2770  NcodeC 12, J2>*1 
2780  CALL  Mmmx<12,  J2> 

2790  IF  Iption*=5  THEN  GOTO  3060 
2800  ! 

2810  !  *****  FIND  CENTER  OF  CIRCLE  ***** 

2820  ! 

2830  Ar < 13, J3)*R3 
2840  Az< 13, J3)*Z3 
2850  NcodeC 13, J3)=l 
2860  CALL  Mn i mx< 13, J3) 

2870  SI  ac*<Z2-Zl  >''<R2-R1  > 

2880  Slbf*-l'Slac 

2890  Slce»<Z3-Z2)/'<R3-R2) 

2900  Sldf*-1/S1ce 
2910  ! 

2920  !  *****  CHECK  FOR  INPUT  ERROR  ***** 

2930  ! 

2940  IF  ABS < S 1 ac -S 1 c e ) > . 00 1  THEN  GOTO  2960 

2950  GOTO  1110 

2960  R4*R1+CR2-R1  / 2 

2970  Z4*Zi*<Z2-Zl)/2 

2980  R5*R2*<R3-R2)/2 

2990  Z5»Z2*<Z3-Z2)/2 

3000  Bbf *Z4-S1 bf *R4 

3010  Bdf *Z5-S1 df *R5 

3020  Rc*<Bbf-Bdf >/<Sldf-S1bf ) 

3030  Zc*S1bf*Rc+Bbf 
3040  Kappa* 1 
3050  GOTO  3090 
3060  Kappa* 2 
3070  Rc  *R3 
3080  Zc*Z3 
3090  Istrt*I 1 
3100  Istp*I2 
3110  J*trt-Jl 
3120  Jstp-J2 
3130  R*trt*Rl 
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Rstp«R2 

Zstrt*Zl 

Zstp-Z2 

CALL  Angl e<Rst rt , Zstrt , Rc , Zc , Angl ) 
CALL  Angl e<Rstp, Zstp, Rc , Zc , Ang2) 

IF  <ABS<Ang2)<=. 00001)  AND  <Angl>PI) 
IF  <ABS<Angl)<=. 00001)  AND  <Ang2>PI) 


THEN 

THEN 


Ang2=2*PI 

Angl*2*PI 


*****  FIND  ANGULAR  INCREMENT  ***** 


Di *ABS< 1st p-Istrt ) 

Dj*ABS< Jstp-Jstrt ) 

I i nc*0 
Ji nc*0 

IF  IstrtOIstp  THEN 
IF  JstrtOJstp  THEN 
Lamdasl 

IF  < I i nc<  >0)  AND  <Jinc<>0) 
Diff«MAX<Di , D j ) 

Iter=Di ff-1 

IF  Lamda=2  THEN  Diff=2*Diff 
IF  Ang2-Angl >PI  THEN  Ang2=Ang2-2*PI 
IF  Ang2-Angl<-PI  THEN  Ang2=Ang2+2*P I 
Delphi =< Ang2-Angl )/Di f f 


Iinc=<Istp-Istrt  >xABS<  Istp-Istrt) 
Jinc=<Jstp-Jstrt)xABS<Jstp-Jstrt) 


THEN  Lamda=2 


*****  CHECK  FOR  INPUT  ERROR  ***** 


IF  <Lamda<>2) 
GOTO  1110 
Io=Istrt 
Jo=Jstrt 


OR  <Di=Dj)  THEN  GOTO  3430 


*****  INTERPOLATE  ***** 


FOR  M= 1  TO  Iter 

IF  Lamda*2  THEN  GOTO  3600 

I=Io+I i nc 

JsJo+J i nc 

CALL  MnimxCI.J) 

NcodeU,  J)*l 

CALL  Ci rc 1 es< Angl , Del  phi , Rstrt, Zstrt,Rc,Zc, I,J) 
IF  Ipt«I3  THEN  GOTO  3770 
Zp<Kp)=Az< I , J) 

Rp<Kp)«Ar< I, J) 

Kp=Kp+l 
GOTO  3770 
I=Io+I i nc 
J*Jo 

Ncode< I , J)*l 
CALL  Mnimx<I,J) 

CALL  Ci rc 1 esCAngl , Del  phi , Rstrt , Zstrt, Rc , Zc , I ,  J) 
IF  Ipt-13  THEN  GOTO  3690 
Zp<Kp)»Az< I , J) 

Rp<Kp)«Ar < I , J) 

Kp*Kp+l 
J* Jo+J i nc 
Ncode< I , J)*l 
CALL  Mnimx<I,J) 

CALL  Ci rc 1 es  < Angl , Del  phi , Rstrt, Zstrt, Rt,Zc, I , J) 
IF  Ipt«13  THEN  GOTO  3770 
Zp<Kp)«Az<I, J) 

Rp<Kp)»Ar< I , J) 

Kp»Kp+l 
Io«I 
Jo»J 
NEXT  M 
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3800  IF  Lamda<>2  THEN  GOTO  3890 
3810  I«Io+Iinc 
3820  Ncode<I,J>«l 
3830  CRLL  Mni mx< I , J) 

3840  CALL  Ci rc 1 es< Angl , De 1 ph i , Rst rt , 2st rt , Rc , 2c , I , J) 

3850  IF  Ipt»13  THEN  GOTO  3890 
3860  2p<Kp)«Az< I , J> 

3870  Rp<Kp>*Ar< I , J> 

3880  Kp-Kp+1 

3890  IF  Kappa»2  THEN  GOTO  1110 

3900  Istrt*I2 

3910  Istp»I3 

3920  Jstrt»J2 

3930  Jstp*J3 

3940  Rstrt*R2 

3950  Rstp“R3 

3960  2strt*22 

3970  2stp=23 

3980  Kappa*2 

3990  IF  Ipt i on=5  THEN  GOTO  4040 

4000  IF  Ipt«13  THEN  GOTO  4040 

4010  Rp<Kp>*R2 

4020  2p<Kp>*22 

4030  Kp=Kp*l 

4040  GOTO  3170 

4050  !  ************************************************************************ 
4060  !  *  CALCULATE  COORDINATES  OF  INTERIOR  POINTS  * 

4070  !  ************************************************************************ 
4080  Inter  ior.'PRINTER  IS  16 

4090  DISP  "  *  Computing  Node  Coordinates  *" 

4100  IF  Max j<=2  THEN  GOTO  4300 

4110  J2=Max j-1 

4120  FOR  N= 1  TO  500 

4130  Resi ds0 

4140  FOR  J»2  TO  J2 

4150  I l=Imin< J)+l 

4160  I2*Imax<J>-l 

4170  FOR  I»I1  TO  12 

4180  IF  Ncode( I , J)-l  THEN  GOTO  4240 

4190  Dr=<Ar<I  +  l,J)+Ar<I-l,J)+Ar<I,J+l)+Ar<I,J-l)>/'4-Ar<I,J> 

4200  Dz*<Az<I+l, J>+Rz(I-l, J)+Az<I, J+l>+Az<I, J-l))/4-Az(I, J> 

4210  Resi d=Resi d+ABS<Dr )+ABS<Dz) 

4220  Ar<I,J)=Ar<I,J)+1.8*Dr 

4230  AzCI, J)=Az<I, J)+1.8*Dz 


4240  NEXT  I 
4250  NEXT  J 

4260  IF  N=1  THEN  Resl*Resid 

4270  IF  ( N*  1 )  AND  <Resid=0)  THEN  GOTO  4300 

4280  IF  Resi d/Resl< . 00001  THEN  GOTO  4300 

4290  NEXT  N 

4300  CALL  Poi nts<#l > 

4310  FOR  Ki«l  TO  25 
4320  FOR  Kj«l  TO  100 

4330  IF  NcodeCKi , K j  >  =  4  THEN  GOTO  4350 

4340  Ncode<Ki , Kj  >=0 

4350  NEXT  Kj 

4360  NEXT  Ki 

4370  FOR  Ki *1  TO  25 

4380  FOR  K  j*l  TO  100 

4390  Ipl»Ki +1 

4400  Jpl»Kj+l 

4410  Ndtest*4 

4420  IF  Ncode<Ki ,  Kj )<  >4  THEN  GOTO  4480 
4430  IF  < Ipl >26)  OR  <Jpl>101>  THEN  GOTO  4480 
4440  IF  Ncode<Ki , Jpl >*4  THEN  Ndtest=2 
4450  IF  Ncode< Ipl , K j  >=4  THEN  Ndtest«l 
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4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 

4770 

4780 

4798 

4808 

4810 

4820 

4830 

4840 

4850 

4860 

4870 

4880 

4890 

4900 

4910 

4920 

4938 

4940 

4950 

4960 

4970 

4980 

4990 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5078 

5080 

5090 

5100 

5110 


IF  (Ncode< Ipl,Kj)=4) 
Ncode<Ki , K j )=Ndtest 
NEXT  Kj 
NEXT  Ki 
1 1 ow* 1 
Ihi gh*25 
J1 ow*l 
Jhi gh*100 
CALL  Plots<#l,F*) 
W-1.5 
BEEP 

WAIT  300/W 
BEEP 

WAIT  150/W 
BEEP 

WAIT  150/W 
BEEP 

WAIT  300/W 
BEEP 

WAIT  600/W 
BEEP 

WAIT  300/W 
BEEP 
PEN 


AND  <Ncode<Ki , Jpl )=4)  THEN  Ndtest=3 


END 
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SUB  Mrii  mx<  I ,  J) 

OPTION  BASE  1 

INTEGER  Ncode, M at erial  ,  Imi  n,  I  m  ax,  J  mi  n,  -J  max,  Code 


Ar<26, 101>,Az<26, 101 >, Ncode <26, 101 > , El em<900, 1 1 ), Mater i al <26, 191 ) 
Nme  1  <  100)  ,  Nemax<  100)  ,  Nemi  n<  100)  ,  Npm  i  r><  100) ,  Npmax<  100) 

Numtc , El emax, Yi el d< 10) 

Imi n< 100), Imax< 100),Jmin<26),Jmax(26), Maxi , Maxj , Nmt 1 , Nbc 
Numnp,Numel , Nummat 

R< 1000), Code< 1000), Xr< 1000), 2< 1000), Xz< 1000) 

Kp,Rp<250),Zp<250) 

Nel e<26) 

1 1 ow, Ihi gh, J 1 ow, Jhi gh 
Npp, Nseg 

Fi le*[15],Line*C81],Que*,P* 

J< Jmi n< I >  THEN  Jmin<I)=J 
THEN  Jmax< I >=J 
THEN  Imin<J)»I 
THEN  Imax<J)=I 


COM 
COM 
COM 
COM 
COM 
COM 
COM 
COM 
COM 
COM 
COM 
IF 

IF  J>Jmax(I) 

IF  I<Imin<J) 

IF  I>Imax<J) 

SUBEXIT 

SUBEND 

SUB  Poi nts<41 ) 

OPTION  BASE  1 

INTEGER  Ncode, Material , Imi n , I  max , Jm i n , Jmax , Code 

COM  Ar<26, 101 >, Az<26, 101 ), Ncode <26, 1 0 1 ) , E 1 em < 900 , 11 ) ,Materi al <26,101) 

COM  Nmel < 100) , Nemax< 100) , Nemi n< 100) , Npmi n< 100) , Npmax< 100) 

COM  Numtc , El emax, Yi e 1 d< 10) 

COM  Imin<100), Imax<100),Jmin<26),Jmax<26),Maxi , Maxj , Nmt 1 , Nbc 
COM  Numnp, Numel , Nummat 

COM  R< 1000) , Code< 1000) , Xr< 1000) ,Z<1000),Xz< 1000) 

COM  Kp, Rp<250) , Zp<250) 

COM  Nel e<26) 

COM  1 1 ow, Ihi gh, J1 ow, Jhi gh 
COM  Npp, Nseg 

COM  Fi le$C15],Line$C81], Que$, PS 
Zero*. 0800001 

!  ************************************************************************ 
!  #  ESTABLISH  NODAL  POINT  INFORMATION  * 

!  ********************** ******************************************** ****** 

Nt  1  *0 

Nodsum*0 

FOR  J»1  TO  Maxj 
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5130  Nstop- 1  max < J) 

5140  FOR  I«Nstart  TO  Nstop 

5150  Nodsum-Nodsum+1 

5160  NEXT  I 

5170  NEXT  J 

5180  N*1sum-0 

5180  J jmax-Max j-1 

5200  FOR  Jj-1  TO  Jjmax 

5210  Nstop=MIN< ImaxC J j ) , Imax<Jj+l))-l 

5220  Nstart=MAX< Imi n< J j  ) , Imin<Jj+l>> 

5230  FOR  li-Nstart  TO  Nstop 
5240  Nel sum*N#1 sum+1 
5250  NEXT  Ii 
5260  NEXT  Jj 
5270  Numnp-Nodsum 
5230  Nuroel *Ne I  sum 
5290  FOR  J-l  TO  Maxj 
5300  Nstart-Imi n< J) 

5310  Nstop-Imax( J) 

5320  FOR  I-Nstart  TO  Nstop 
5330  Np*FNNod*< I , J) 

5340  R  (  Np  >  -Ar  <  I ,  J  > 

5350  2<Np)=Az< I ,  J) 

5360  NEXT  I 
5370  NEXT  J 

5380  !  ************************************************************************ 
5390  !  *  READ  AND  ASSIGN  BOUNDARY  CONDITIONS  * 

5400  !  *********»**************)HHHHHHHf*********JHHHHHf******#*»»**iHf»****jHf**** 

5410  !  *  INITIALIZE 

* 

5420  !  ************************************************************************ 
5430  FOR  1-1  TO  Numnp 
5440  Cod»<I>*0 

5450  IF  ABS<<R< I ><=Zero)  AND  <Npp=0>>  THEN  Code(I)=l 
5460  Xr(l)»0 

5470  Xz<I)-0 

5480  NEXT  I 

5490  IF  Nbc *0  THEN  GOTO  5680 
5500  FOR  Ibcon-1  TO  Nbc 
5510  READ  #1 ; L i ne* 

5520  Il-VAL(Lin*»C2,6]) 

5530  I2-VAL<Un«$[7,  11]) 

5540  Ji-VAL<Lin»$tl2, 16]) 

5550  J2-VAL<Lin«*C17,21  ]) 

5560  !  Con«VAL<Lin**C22,31]) 

5570  !  Rcon-VAL<Lin«*C32,41]) 

5580  <  Zcon«VAL<Lin«*C42,51 ]> 

5590  FOR  I-Il  TO  12 
5600  FOR  J-Jl  TO  J2 
5610  Np-FNNode< I , J) 

5620  Cod*<Np)*Con 
5630  Xr<Np)«Rcon 
5640  Xz<Np)»Zcon 
5650  NEXT  J 
5660  NEXT  I 
5670  NEXT  Ibcon 
5680  FOR  J-1  TO  Maxj 
5690  Nstart-Imi n<J) 

5700  Nstop-Imax< J) 

5710  FOR  I -Nst art  TO  Nstop 
5720  Np-FNNod*< I ,  J) 

5730  Ncod«< I , J)-4 
5740  Ar < I , J)*R<Np> 

5750  Az<I, J)«Z<Np) 

5760  NEXT  I 
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5770  NEXT  J 
5780  !  Nbc»Nbcl 
5790  SUBEXIT 
5800  SUBEND 

5810  SUB  Angl e<R, Z, Rc , Zc , Ang) 

5820  OPTION  BASE  1 

5830  INTEGER  Ncode, Mater i al , Imi n, Imax, Jmi n, Jmax, Code 

5840  COM  Ar<26, 101 ), Az<26, 101 ), Ncode <26, 10 1 > , E 1 em<900, 1 1 >, Material <26, 101 > 

5850  COM  Nme  1  < 100) , Nemax< 100 ) , Nemi n< 1 00) , Npmi n< 100) , Npmax < 100) 

5860  COM  Numt  c , E 1 emax, Y i e 1 d< 1 0) 

5870  COM  Im i n< 100) , I  max  < 100),Jmin<26>, Jraax  <  26 ) , Max i , Max j , Nmt 1 , Nbc 

5880  COM  Numnp, Numel , Nummat 

5890  COM  R(1000),Code<1000),Xr<1000),Z<1000),Xz(1000) 

5900  COM  Kp,Rp<250),2p<250) 

5910  COM  Nel e<26) 

5920  COM  1 1 ow, Ihi gh, J1 ow, Jhi gh 

5930  COM  Npp, Nseg 

5940  COM  Fi lePC15I, Line $1813, QueP,PP 

5950  Dl*Z-Zc 

5960  D2»R-Rc 

5970  IF  flBS<R-Rc)>. 00000001  THEN  GOTO  6050 

5980  Ang*PI/2 

5990  IF  Dl>. 00000001  THEN  SUBEXIT 

6000  Ang*-Ang 

6010  SUBEXIT 

6020  ! 

6030  !  *****  ALLOW  CIRCLE  TO  CROSS  AXIS  ***** 

6040  ! 

6050  Ang*ATN<DlxD2> 

6060  IF  <D1>=0)  AND  <D2<0)  THEN  Ang=Ang+PI 

6070  IF  <D1<=0)  AND  <D2<0)  THEN  Ang=Ang-PI 

6080  SUBEXIT 

6090  SUBEND 

6100  DEF  FNNode< I , J )  «  *****  THIS  IS  THE  SUBFUNCTION  Node  ***** 

6110  OPTION  BASE  1 

6120  INTEGER  Ncode, Mat eri al , Im i n, Imax, Jmi n, Jmax, Code 

6130  COM  Ar<26, 101),Az<26, 101 ) , Ncode <26, 101 ) , El em<900, 11 ) , Mat eri al <26,101) 
6140  COM  Nmel < 100) , Nemax< 100) , Nerai n< 100) , Npmi n< 100) , Npmax< 100) 

6150  COM  Numtc,Elemax,Yield<10) 

6160  COM  Imi n< 100) , Imax < 100),Jmin<26),Jmax<26>, Maxi , Max j , Nmt 1 , Nbc 
6170  COM  Numnp, Nume 1 , Nummat 

6180  COM  R<1000),Code<1000),Xr<1000),Za000),Xz<1000) 

6190  COM  Kp , Rp  <  250 ), Zp  <  250 ) 

6200  COM  Nel e<26) 

6210  COM  1 1 ow, Ihi gh, J1 ow, Jhi gh 

6220  COM  Npp, Nseg 

6230  COM  Fi l eSI 151 , L i net t 813 , Que* ,P$ 

6240  Node«0 

6250  FOR  J j » 1  TO  J 

6260  Nstart*Imi n< Jj ) 

6270  Nstop«Imax< Jj ) 

6280  FOR  Ii^Nstart  TO  Nstop 

6290  Node-Node+1 

6300  IF  <Jj*J>  AND  < I i = I )  THEN  RETURN  Node 

6310  NEXT  li 

6320  NEXT  Jj 

6330  RETURN  Node 

6340  FNEND 

6350  SUB  P1ots<«l,F$) 

6360  OPTION  BASE  1 

6370  INTEGER  Ncode, Mater i al , Imi n, Imax, Jmi n, Jmax, Code 
6380  PLOTTER  IS  "GRAPHICS" 

6398  GRAPHICS 
6400  BEEP 

6410  IF  P»«"Z"  THEN  GOTO  6530 

6420  INPUT  "Do  you  want  to  create  a  Grid  Fi 1 e?< Y/N) " , RP 


6430  IF  R$»"N"  THEN  GOTO  6460 
6440  DIM  Gr  i d<2000, 3) 

6430  PLOTTER  IS  Grid<*) 

6460  INPUT  "Hard  copy <Y/N) " , flnsl 
6470  IF  Ans$*“N"  THEN  GOTO  6530 
6480  PEN  1 

6490  OUTPUT  7,5;  MIP1000,0,  15102,  11400'* 

6500  PLOTTER  IS  "9872A" 

6510  PRINT  "Prepare  platten  for  LARGE  paper  then  C  “ 

6511  PRINT  11  0  " 

6512  PRINT  “  N  " 

6513  PRINT  "  T  " 

6520  PAUSE 

6521  PRINT  PAGE 

6530  COM  Ar<26, 101 ) , Az<26, 101 ) , Nc ode <26, 101 ) , E 1 em<9O0, 1 1 >, Materi al <26,101) 
6540  COM  Nme 1 < 1 00 1 , Nemax < 1 00 > , Nero i n < 1 00 > , Npm i n < 1 00 ) , Npmax < 1 00 ) 

6530  COM  Numt  c  ,  E 1  e.^ax ,  Y  i  e  1  d <  1 0) 

6560  COM  Imi n( 100) , ifnax<100),Jmin<26),Jmax<26),Maxi , Max j , Nmt 1 , Nbc 
6570  COM  Numnp , Nume 1 , Nummat 

6580  COM  R< 1 000) , Code< 1 000) ,Xr<1000),2<1000),Xz(1000) 

6590  COM  Kp,Rp<250),2p<250) 

6600  COM  Ne 1 e  <  26 ) 

6610  COM  1 1 ow, Ihi gh, J1 ow, Jhi gh 
6620  COM  Npp,Nseg 

6630  COM  File*n5],L»ne$C813,Que*,Pt 

6640  Scale:!  ******  AUTOMATIC  SCALING  ****** 

6650  Ihf  1 g*l 

6660  1 1 f 1 g*l 

6670  FOR  k*Ilow  TO  Ihigh 

6680  IF  Ncode<K, Jhi gh)=0  THEN  GOTO  6710 

6690  2max*MAX< 2max , Az < K , Jh i gh ) ) 

6700  Ihf 1 g*0 

6710  IF  Ncode<K, Jlow)=0  THEN  GOTO  6740 
6720  2min=MIN<2min,Az<K,Jlow)) 

6730  1 1 f 1 g*0 

6740  NEXT  K 

6750  IF  Ihf 1 g*l  THEN  GOTO  6930 

6760  IF  Ilflg-1  THEN  GOTO  6950 

6770  Ihflg*l 

6780  1 1 f 1 g*l 

6790  FOR  K*Jlow  TO  Jhigh 

6800  IF  Ncode< Ihi gh, K)=0  THEN  GOTO  6830 

6810  Rmax=MAX<Rmax, Ar< Ihi gh, K) ) 

6820  Ihflg-0 

6830  IF  Ncode< 1 1 ow, K)=0  THEN  GOTO  6860 
6840  Rmi n=MIN<Rmi n, Ar < 1 1 ow, K> ) 

6850  II f 1 g*0 
6860  NEXT  K 

6870  IF  Ihf 1 g*l  THEN  GOTO  6970 
6880  IF  1 1 f 1 g* 1  THEN  GOTO  6990 
6890  CALL  Full <2max,A,B,D) 

6900  IF  2max<  =Rmax/ .  9  THEN  CALL  Rscal  e<Rw>ax,  A,  B,  D) 

6910  SCALE  A, 6, 0,11 
6920  GOTO  7010 
6930  Jhigh*Jhigh-l 
6940  GOTO  6650 
6950  J 1 ow* J 1 ow+ 1 
6960  GOTO  6650 
6970  Ihi gh*Ihi gh- 1 
6980  GOTO  6770 
6990  I1ow*Ilow+l 
7000  GOTO  6770 

7010  !  ******  PLOT  IN  THE  J- DIRECTION  ****** 

7020  Movflg*! 

7030  Jhml*Jhigh 
7040  INTEGER  Ki ,Kj 


39 


7050 

7060 

7070 

7080 

7090 

7100 

7110 

7120 

7130 

7140 

7150 

7160 

7170 

7180 

7190 

7200 

7210 

7220 

7230 

7248 

7250 

7260 

7270 

7280 

7290 

7300 

7310 

7320 

7330 

7340 

7350 

7360 

7370 

7380 

7390 

7400 

7410 

7420 

7430 

7440 

7450 

7460 

7470 

7480 

7490 

7500 

7510 

7520 

7530 

7540 

7550 

7560 

7570 

7580 

7590 

7600 

7610 

7620 

7630 

7640 

7650 

7660 

7670 

7680 

7690 

7700 


FOR  Ki»llou  TO  Ihigh 
FOR  Kj*Jl ow  TO  Jhml 
Ndcode»Ncode<Ki , K j ) 


DIRECTION  ****** 


IF  (Ndcodes3)  OR  <Ndcode=2>  THEN  GOTO  7110 
Movf  1 g=l 
GOTO  7230 

IF  Movf 1 g» 1  THEN  GOTO  7160 
Xzj«Az<Ki , Kj  +  1 > 

Yr i ■RrCKi , Kj+1 > 

DRAW  Xzj , Yr i 
GOTO  7230 
Movf 1 g»0 
Xzj«Az<Ki , Kj ) 

Yri=Ar<Ki ,Kj> 

MOVE  Xzj , Yr i 
Xz j*Az<Ki , K j+1 ) 

Yr i =Ar <Ki ,Kj+l> 

DRAW  Xzj , Yr i 
NEXT  Kj 
Movf lg*l 
NEXT  Ki 

!  ******  PLOT  IN  THE  I-DIRECTION  ****** 

FOR  Kj*Jlow  TO  Jhigh 
Ihml=Ihi gh 
FOR  Ki=I1ow  TO  Ihml 
Ndcode-NcodeCKi ,Kj) 

IF  <Ndcode=3>  OR  <Ndcode=l )  THEN  GOTO  7340 
Movf 1 g«l 
GOTO  7460 

IF  Movf 1 g« 1  THEN  GOTO  7390 
Xzj*Az<Ki+l,Kj) 

Yri *Ar<Ki +1 , Kj  > 

DRAW  Xzj, Yri 
GOTO  7460 
Movf 1 g«0 
Xzj=Az<Ki , K j ) 

Yri *Ar<Ki ,Kj) 

MOVE  Xzj, Yri 
Xzj*Az<Ki +1 , Kj  ) 

Yri =Rr <Ki +1 , Kj ) 

DRAW  Xzj, Yri 
NEXT  Ki 
Movf 1 g«l 
NEXT  Kj 

IF  Ans**'^"  THEN  OUTPUT  7,5;"SP0" 

IF  An»$B"Y"  THEN  PLOTTER  7,5  IS  OFF 
IF  R$«"N"  THEN  GOTO  7580 
PLOTTER  Gr i d<*>  IS  OFF 

INPUT  "Enter  name  of  grid  file:  G _ ",Gfi 

Gf  i  le»*Gf  i  le*8,":C12" 

CREATE  Gf i le«,300 
ASSIGN  *2  TO  Gf i le* 

PRINT  *2; Gr i d(«) 

DIM  El ement < 100, 3) 

INTEGER  M*teriol 

DIM  Met er i o 1 <M*x i , Max j ) 

PRINT  PAGE 

PRINT  "When  you  are  finished  vie 


,  Gf  i  1  e$ 


PRINT  "When 
PRINT  " 
PRINT  " 
PRINT  " 
PAUSE 
GCLEAR 
PRINT  PAGE 
IF  Que*«"Y" 


are  finished  viewinc 


THEN  PRINTER  IS  0 


FOR  Mat< 


TO  Nmt 1 


7710 

7720 

7730 

7740 

7750 

7760 

7770 

7780 

7790 

7800 

7810 

7820 

7830 

7840 

7850 

7860 

7870 

7880 

7890 

7900 

7910 

7920 

7930 

7940 

7950 

7960 

8030 

8040 

8050 

8060 

8070 

8080 

8090 

8100 

8110 

8120 

8130 

8140 

8150 

8160 

8170 

8180 

8190 

8260 

8210 

8220 

8230 

8240 

8250 

8260 

8270 

8280 

8290 

8300 

8310 

8320 

8330 

8340 

8350 

8360 

8370 

8380 

8390 

8400 

8410 

8420 


READ  4l;Line* 

Hue*VAL<Line*Cl,5]> 

I l*VAL<Line*[6,  10]) 

Ih*VAL<Line$Cll, 151) 

J1  =VAL<Li ne$C  16, 20  3  > 

Jh=VAL<L ine$[21,253> 

IF  Que**"Y“  THEN  PRINT  USING  I  mage 3; Hue, II , Ih,  J1 ,  Jh 
Image3:  IMAGE  5<5D> 

FOR  J=J1  TO  Jh-1 

FOR  1*11  TO  Ih-1 

IF  R  ♦  = "  S "  THEN  GOTO  7830 

IF  P*="A"  THEN  GOTO  7830 

Materi al  <  I , J)*Hue 

Mat er i ol <  I ,  J  > =Hue 

PLOTTER  IS  E 1 ement  <  *  > 

MOVE  Az< I , J) , Ar< I , J) 

PEN  1 

DRAW  Az<I+l,J),Ar<I+l,J> 

DRAW  Az<I+l, J+l),Ar<I+l, J+l) 

DRAW  AzCI, J+l),Ar<I, J  +  l) 

DRAW  AzCI, J),Ar<I, J) 

PLOTTER  E 1 ement  <  * )  IS  OFF 
Huey*Hue 

IF  Hue >6  THEN  Huey=Hue-4.5 

AREA  COLOR  Hueyx6, 1 , 1 

MAT  PLOT  Element, FILL 

NEXT  I 

NEXT  J 

NEXT  Mat 

PRINTER  IS  16 

INPUT  "Do  you  want  to  create  a  Supergrid  Fi 1 e?< YxN) " , Rf 
IF  R  *  * " Y "  THEN  GOTO  Supergrid 
GOTO  8560 

Supergr i d: ! 

I 

!  *********  READ  MATERIAL  INFORMATION  CARDS  ********** 

j 

READ  4i;Line* 

FOR  1*1  TO  Nummat 
READ  4l;Line* 

Numid=VAL<Line*Cl,53) 

READ  4l;Line* 

READ  41 ; Li ne* 

Yiel d<Numi d) =VAL < L i ne$C 3 1 , 40]  > 

NEXT  I 

FOR  J*1  TO  Max j 
FOR  1=1  TO  Maxi 

IF  MateriaKI,  J)=0  THEN  GOTO  8460 

j 

!  ****  Exchange  mat'l  4  for  Yield  Strength  **** 

j 

FOR  L*1  TO  10 

IF  Mater i al < I , J)=L  THEN  Mater i al (I , J)=Yi el d<L> 

NEXT  L 

! 

!  ****  Create  Element  Information  Array  **** 

i 

E 1 num=E 1 num+ 1 

Elem(E1num,l)=Elnum! . . . Element  4 

E 1  em<E  1  num,  2>=Mat  er  i  al  <  I ,  J )  ! . Material  Yield  Strength 

E 1  em<E  1  num,  3 )=Az<  I ,  J )  ! . LL-z 

El  em<El  num,  4 )*Ar  <  I ,  J )  ! . LL-r 

El em<El num,  5)=Az<  1  + 1 ,  J)  > . UL-z 

El  em<E  1  num,  6)=Ar <  I*  1 ,  J)  ! . UL-r  UL _ UR 

El  em(El  num,  7>*Az<  1  +  1 ,  J+l  >  ! . .UR-z  | _ |< Element 

E 1  em<E  1  num,  8) *Ar  <  I  + 1 ,  J+ 1 ) ! . UR-r  LL  LR 


8430 

8440 

8450 

6460 

8470 

8480 

6490 

8500 

8510 

8520 

8530 

8540 

8550 

8560 

8570 

8580 

8590 

8600 

8610 

8620 

8630 

8640  Imagel 

8650 

8660 

8670 

8680 

OTO  8710 

8690 

8700 

8710  Image2 
8720 
8730 
8740 
8750 
8760 
8770 
8780 
8790 
8800 
8810 
8820 
8830 
8840 
8850 
8860 
8870 
8880 
8890 
8900 
8910 
8920 
8930 
8940 
8950 
8960 
8970 
8980 
8990 
9000 
9010 
9020 
9030 
9040 
9050 
9060 
9070 


El  em<E1  num,  9)=Az<I,J+l)!. . LR-z 

El  em(E1  num,  10)*Ar<  I ,  J+l  >  !  . .  . . LR-r 

El  em(El  num,  ll)=Materio1<I,J)! . Material  I.D.# 

NEXT  I 
NEXT  J 

El emax=E 1 num 

IF  R*="N"  THEN  GOTO  8580 

INPUT  "Enter  name  of  Supergrid  file:  S _ ",Sfile* 

Sf ile*-Sf i  le$&":C12" 

CREATE  Sf i 1 e* , E 1 emax* 1 1 /4+ 1 
ASSIGN  #2  TO  Sf i 1 e$ 

PRINT  #2 ; E 1 emax 
MAT  PRINT  42; El em 
BEEP 

PRINTER  IS  16 

PRINT  "Would  you  like  a  nodal  point  1 i st i ng?< YxN) " 

INPUT  A$ 

IF  A*="N"  THEN  GOTO  8780 

INPUT  "Does  the  grid  contain  the  I,J  or i g i n? < Yxn ) " , 

PRINTER  IS  0 
PRINT  USING  Image  1 

:  IMAGE  3X"r,,3X"J"  ,7X"R",  10X"2",  10X"NP" 

FOR  K= 1  TO  Max j 
FOR  L*1  TO  Maxi 

IF  <D$="N">  AND  <  K=  1  >  AND  <L=1)  THEN  GOTO  8710 

IF  ((KOI)  OR  <  L<  >  1 )  >  AND  <Ar<L,K)=0)  AND  <Az<L,K)=0)  THEN  G 

J i m= J i m+ 1 

PRINT  USING  Image2;L,K,Ar(L,K),Az(L,K),Jim 
:  IMAGE  1X,DDD, 1 X , DDD, 5X, DD . DDD, 5X , DD . DDD , 5X , DDDD 

NEXT  L 
NEXT  K 
BEEP 

WAIT  200 
BEEP 

WAIT  200 
BEEP 
SUBEXIT 
SUBEND 

SUB  Ci rc 1 esKAngl , Del  phi ,Rstrt,Zstrt,Rc,Zc, I , J> 

OPTION  BASE  1 

INTEGER  Ncode,Material , Imin, I  max , Jm i n, Jmax , Code 

COM  Ar<26, 101),Az<26, 1 0 1 ) , Nc ode < 26 , 1 0 1 ) , E 1 em < 900 , 11), Material <26,101) 
COM  Nme l<100),Nemax<100),Nemin(100),Npmin(100), Npmax <100) 

COM  Numtc , El emax, Yi e 1 d< 10) 

COM  Imin(100), Imax<100),Jmin<26),Jmax<26),Maxi ,Maxj,Nmt 1 ,Nbc 
COM  Numnp, Nume 1 , Nummat 

COM  R< 1000) , Code < 1000) , Xr< 1000), Z< 1000) , Xz< 1000) 

COM  Kp,Rp<250),Zp(250) 

COM  Nel e<26) 

COM  1 1 ow, Ihi gh, J1 ow, Jhi gh 
COM  Npp,Nseg 

COM  Fi le*C15],Line*[813,Que*,P$ 

!  ************************************************  ! 

!  FIND  INTERSECTION  OF  LINE  8,  CIRCLE  =  NEW  R  &  Z 
!  ************************************************  ! 

Angl«Angl+Delphi 

Rr*SQR<  <Rstrt-Rc  )-'2+<Zstrt-Zc  )-2) 

Ar < I , J)»Rc+Rr*COS<Angl ) 

Az< I , J)»Zc+Rr*SIN(Angl ) 

SUBEXIT 

SUBEND 

SUB  Ful 1 (Zmax, A, B, D) 

A*Zmi n 
B*Zmax 
D*Zmaxxl . 23 


|  9080 

SU8EXIT 

1  9090 

SUBEND 

I  9100 

SUB  Rscal e<Rmax, ft,  B,  D> 

R  9110 

A*Zmi  n 

I  9120 

B*Rmax*l . 23*1 . 2 

|  9130 

D*Rmax*l .  2 

g  9140 

SUBEXIT 

I  9150 

SUBEND 

1  9160 

SUB  Gridl<#2,Gfi let, Grid) 

£  9170 

PLOTTER  IS  "GRAPHICS" 

1  9180 

GRAPHICS 

j  9190 

SUBEXIT 

1  9200 

SUBEND 

|  9210 

SUB  Gr i d2<#2, Gf i 1 e$, Gr i d) 

1  9220 

INPUT  "Enter  name  of  grid  f 

t  9230 

Gfi  1 e$*Gf i  1 e$&" : C12" 

i  9240 

CREATE  Gfi le*,400 

[  9250 

ASSIGN  «2  TO  Gfi 1e* 

>  9260 

DIM  Gr i d( 1500, 3> 

S  9270 

PRINT  #2} Gr i d(*> 

I  9230 

SUBEXIT 

r  9290 

■41 

SUBEND 

1  !  *********  PROGRAM  " PREP1 "  OF  THE  PREPPY  PACKAGE  ********** 

10  OPTION  BASE  1 

20  PRINTER  IS  16 

30  PLOTTER  IS  "GRAPHICS" 

40  GRAPHICS 

50  CALL  Plots 

100  END 

110  SUB  Plots 

120  OPTION  BASE  1 

130  INPUT  “Enter  name  of  Supergrid  f i 1 e : " , Sf i 1 s$ 

140  DISP  “*  Reading  in  Supergrid  File:  “ ;Sfile* 

150  Sf  i  1e»=Sf  i  1  e*8." :  C12" 

160  ASSIGN  ttl  TO  Sf i 1 e* 

170  DIM  El em<850, 11) 

180  READ  ttl ; E 1 emax 

190  REDIM  El em<El emax, 1 1 ) 


200 

210 

211 

220 

221 

230 

240 

250 

251 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 


MAT  READ  »l;E1em 
MAT  SEARCH  E 1 em < * , 7 > , MAX ; Zmax 1 
MAT  SEARCH  E 1 em < * , 3 > , M I N ; 2m i n 1 
MAT  SEARCH  E 1 em < * , 9 > , MAX ; Zmax2 
MAT  SEARCH  E 1 em ( * , 5 > , M I N ; Zm i n2 
MAT  SEARCH  E 1 em < * , 8 > , MAX ; Rmax 1 
MAT  SEARCH  E 1 em < * , 1 0 ) , MAX ; Rmax2 
Zmax  =  MAX ( Zmax 1 ,  Zmax2> 

Zmi n=MIN<Zmi nl , Zmi n2> 


Rmax=MAX ( Rmax 1 , Rmax2) 

DIM  Stress$<850)C803 
DIM  El ement < 10, 3) 

SHORT  Stress<850,7) 

REDIM  Stress<El emax, 7) 

INPUT  "Material  plot  or  Stress  p 1 ot ? (MxS) " , Choi  cel 
IF  Choi ce*="M"  THEN  GOTO  Scalemenu 
INPUT  "Enter  name  of  output  file:" ,Out* 

Out*  =  Out*8,":C12" 

ASSIGN  tt3  TO  Out ♦ 

DISP  "  *  Organizing  Stress  Array 

READ  #3; Stress$< 1 ) 

FOR  1=1  TO  Elemax 
READ  tt3 ; St resstt ( I ) 

Stress< I , 1 >=VAL<Stress#< I > [ 1 , 43  > 
Stress<I,2)=VAL<Stress*<I)C5, 111) 
Stress<I,3)=VAL<Stress*<I)C13,  19  3) 

St ress< I, 4 )=VAL <Stress*<I >121,273) 

Stress<I, 5) =VAL (Stress*  <I)C29,353> 
Stress<I,6>=VALCStress$<I>C37,433> 
Stress<I,7>=Elem<I, 11) 


NEXT 

Scalemenu:  ! 


I 


490 

PRINT  "Scale  Menu: 

FULL 

500 

PRINT  " 

1x2 

510 

PRINT  " 

2x2 

520 

PRINT  " 

1X4 

530 

PRINT  " 

2x4 

540 

PRINT  " 

3x4 

550 

PRINT  " 

4x4 

560 

INPUT  Menu* 

570 

PRINT  PAGE 

580 

GCLEAR 

590 

IF  <Zmax>Rmax) 

AND 

<Menu$=" 

600 

IF  <Zmax>Rmax) 

AND 

<Menu*=" 

610 

IF  <Zmax>Rmax> 

AND 

<Menu*=" 

620 

IF  (Zmax)Rmax) 

AND 

<Menu$=" 

630 

IF  (Z»ax>Rmax> 

AND 

<Menu*=" 

640 

IF  <Zmax>Rmax) 

AND 

<Menu»=" 

650 

IF  <Zmax>Rmax) 

AND 

<Menuf=" 

660 

IF  . 99#Zmax< *Rmax 

THEN  CALL 

FULL")  THEN  CALL  Fu 1 1 < Zmax , A , B , D ) 

1^2")  THEN  CALL  One_hal f (Zmax, A, B, D) 
2^2")  THEN  CALL  Two_hal f ( Zmax , A , B , D ) 
1x4")  THEN  CALL  One_f our t h ( Zmax , A , B , D ) 
2x4")  THEN  CALL  Two_f ourt h < Zmax , A , B , D ) 
3x4")  THEN  CALL  Thr ee_f ourt h < Zmax , A , B , D ) 
4x4")  THEN  CALL  Four_f  our  t.  h  <  Zmax  ,  A ,  B ,  D  ) 
Rscale<Rmax,A,B,D> 


661 

6?0 

680 

690 

700 

701 
710 
720 

D  (Di r#<  >" 
730 

D  (Dir$<>" 
740 

D  ( D i r$<  > " 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850  Effect 

851 
860 
870 
880 
900 
910 
915 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 

1  160  ! 

1170  ! 

1180  ! 

1190 
1200 
1210 
1220 
1230 
1240 
1250  ! 

1260  ! 


SCALE  A, 6, 0, D 

IF  Choi ce$="M"  THEN  GOTO  Effective 

PRINT  “Specify  Stress  Component:  R ,  2 .  T ,  RZ ,  EFF  11 

INPUT  Dir$ 

IF  (Menu$="FULL" )  OR  (Menu*=“ 1/2" )  OR  ( Menu*® " 1 /4 " )  THEN  A=Zmin 
IF  Chart f 1 ag  = 1  THEN  SCALE  A,B,0,D 

IF  (Dir*<>"R")  AND  (Dir*<>"Z")  AND  (Dir*<>"T")  AND  (Di r$< >"RZ" ) 
")  THEN  PRINT  "INVALID  ENTRY" 

IF  (Di r*<  >"R" )  AND  <Dir$<>"Z">  AND  (Dir*<>"T")  AND  (Dir*<>"RZ") 
">  THEN  OUTPUT  7 , 6 ; " BP7 , 300 , 5 ; BP6 , 300 , 5 " 

IF  <Dir*<>"R“>  AND  (Dir*<>"Z")  AND  (Di r*< >"T" )  AND  (Di r$<>"RZ" ) 

a  \  riiru  ATrt  y-rtrt 


i  r  i  ncn  ountc  n  , 

IF  (Di  r*<  >"R"  )  AND  (Dir*<>"Z“)  AND  (Dir*<>"T") 
•")  THEN  PRINT  "INVALID  ENTRY" 

IF  (Di r*< >" R" >  AND  (Dir»<>"2">  AND  (Di r$< >"T" > 
THEN  OUTPUT  7 , 6 ;  " BP7 , 300 , 5 ; BP6 , 300 , 5 " 

IF  (Di r$<  >"R" )  AND  (Dir$<>"Z">  AND  (Di r$< >"T" > 
•“>  THEN  GOTO  690 


EFF")  THEN  OUTPUT  7 , 6 ; " BP7 , 300 , 5 ; BP6 , 300 , 5 " 

IF  (Di r$<  >"R" >  AND  (Dir$<>"Z")  AND  (Di r$< >"T" )  AND  (Dir$<>"RZ") 
EFF")  THEN  GOTO  690 

IF  Dir*<>"EFF"  THEN  GOTO  790 
PRINT  PAGE 

PRINT  "Component  =  EFF" 

GOTO  Effective 

PRINT  "Specify  Di  r  ect  i  on:  (-/+) 11 
INPUT  Ct$ 

PRINT  PAGE 
GCLEAR 

PRINT  " Component  =  “ ; D i r$, " D i rec  t i on  =  ";Ct$ 

GOTO  Other_di rect i on 
ive:  IF  Choi ce$="S"  THEN  CALL  Chart 

IF  (Menut="FULL" )  OR  ( Menu*= " 1 /2 " )  OR  ( Menu*= " 1 / 4 " )  THEN  A=Zm 

SCALE  A , B , 0 , D 

MAT  SORT  Stress(*,l) 

FOR  1=1  TO  Elemax 
Hue=El em(  1 ,  11) 

IF  Hue=l  THEN  1240!***IF  YOU  WANT  AIR  PLOTTED  COMMENT  THIS  LI 
IF  Cho i c e$= " M "  THEN  GOTO  960 
Intensity=Stre£s(I,6)/Elem(I,2) 

Intense=Intensi ty 
IF  Intensity>l  THEN  Intensity=l 
IF  lntensity<0  THEN  Intensity®© 

IF  Hue >6  THEN  Hue  =  Hue~4.5 
Hue=Hue/6 

PLOTTER  IS  E 1 ement ( * ) 

MOVE  Elem(I,3),Elem(I,4) 

PEN  -1 

DRAW  El em< I , 5) , El em< I , 6) 

DRAW  El em( I , 7) , El em( I ,8) 

DRAW  Elem(I,9),Elem(I, 10) 

DRAW  Elem(I,3),Elem(I,4) 

PLOTTER  E 1 ement ( * )  IS  OFF 
IF  Choice$="M"  THEN  Intensity=l 
AREA  COLOR  Hue, Intensi ty,  1 
MAT  PLOT  Element, FILL 
IF  Choi ce$="S"  THEN  1170 
MOVE  Elem(I,3),Elem(I,4) 

PEN  -1 

DRAW  Elem(I,5),Elem(I,6) 

DRAW  Elem(I,7),Elem(I,8) 

DRAW  El em( I , 9) , El em( I , 10) 

DRAW  Elem(I,3),Elem(I,4) 


*****  ELEMENTS  ABOVE  YIELD  ***** 


IF  ( I nt ense< 1 . 02 )  OR  <Dir*< 
MOVE  Elem( I , 3) , E 1 em( I, 4) 
DRAW  El em( I , 7) , El em( I ,8) 
MOVE  Elem(I,5),Elem(I,6) 
DRAW  Elem(I,9),Elem(I, 10) 
NEXT  I 


******************************** 


'EFF")  THEN  GOTO  1240 


1270  ! 

1280  PRINTER  IS  16 

1290  IF  Choi ce$="M"  THEN  INPUT  "Look  at  another  sect i on?<Y/N> " , Lo 

ok* 


1300  IF  Look$="N"  THEN  3270 

1310  IF  Choi ce*="M"  THEN  GOTO  Scalemenu 

1320  INPUT  "More  data  ( Y/N> “ , Huh* 

1330  IF  Huh*«"N"  THEN  GOTO  3270 

1340  PRINT  PAGE 

1350  PRINT  "Stress  Component  or  Output  File:  SCxQP" 

1360  INPUT  An* 

1370  IF  <An*<  >"SC" )  AND  <An$O"0P">  THEN  PRINT  "INVALID  ENTRY" 

1380  IF  <An*<>"SC">  AND  <An*O”0P")  THEN  GOTO  1350 

1390  PRINT  PAGE 

1400  IF  An*="SC"  THEN  GOTO  Scalemenu 

1410  IF  An*="OP"  THEN  GOTO  330 

1420  Other  direction:! 

1440  DIM  One<350,5>,Two<850,5),Three<300,5>,Four<850,5),Fi ve<300,5>,Six<850, 
5)  ,  Seven <850, 5)  ,  Ei  ght  <  300, 5), Nine  <300, 5), T en<  850, 5>,Maxva1s<10,8> 

1450  Kl=K2=K3*K4=K5=K6=t<7=K8=K9=K10=0 

1460  MAT  SORT  Stress<*,7>  DES 

1470  Loc=Elemax+l 

1480  FOR  Mat  1 *1  TO  10 

1490  MAT  SEARCH  St ress < * , 7 > , LOC < =Mat 1 > ; 0 1 d 1 oc 

1500  IF  01 dl oc=El emax+1  THEN  1520 

1510  ON  Stress<Loc-l , 7)  GOTO  One, Two , Three , Four , F i ve, Si x , Seven, Ei ght , Ni ne, Te 

n 


1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

I860 

1870 

1880 

1890 

1900 


!  CONTINUE 
NEXT  Mat  1 
GOTO  3140 

One: !  ************************************** 

REDIM  0ne<Loc-01 dl oc , 5) 

FOR  1=01 dl oc  TO  Loc-1 
K 1 =K 1 + 1 
FOR  J*1  TO  5 
One<Kl ,J)=Stress<I,J> 

NEXT  J 
NEXT  I 
Loc=01 dl oc 
FOR  K=1  TO  4 

MAT  SEARCH  One < *, K  + 1 ) , MAX; Maxval s < 1 , K ) 

MAT  SEARCH  One<», K+l ), MIN; Maxval s< 1 , K+4> 

NEXT  K 
Maxval mark=l 

CALL  Di  rec  t i onpl ot<Elem<*>, One  <*),Kl,Dir$,Ct$, Maxval mark , Maxval s< * > ) 
GOTO  1520 

Two:  !  ************************************ 

REDIM  Two<Loc-01 dl oc , 5) 

FOR  1=01 dl oc  TO  Loc-1 

K2-K2+1 

FOR  J*1  TO  5 

Two<K2,J)=Stress<I,J) 

NEXT  J 
NEXT  I 
Loc*01 dl oc 
FOR  K* 1  TO  4 

MAT  SEARCH  Two<*, K+l ), MAX; Maxval s<2, K) 

MAT  SEARCH  Two<*, K  + 1 >, MIN; Maxval s <2, K+4 ) 

NEXT  K 

Maxval mark  =  2 

CALL  Directionplot<Elem<*),Two<*),K2,Dirf,Ctf,Maxvalmark,Maxva1s<*)> 

GOTO  1520 

Three: !  ************************************ 

REDIM  Three <Loc -01 dl oc , 5> 

FOR  1=01 dl oc  TO  Loc-1 
K3-K3+1 
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1918 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

) 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2168 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2548 

2550 


FOR  J»1  TO  5 

Three <K3, J)*Stress<I,J) 

NEXT  J 
NEXT  I 
Loc*01 dl oc 
FOR  K«1  TO  4 

MAT  SEARCH  Three<*, K+l > , MAX; Maxval s<3, K> 

MAT  SEARCH  Three<*f K+l ) , MIN; Maxval s<3, K+4) 

NEXT  K 

Maxval mark=3 

CALL  Directi onpl ot<Elem<*),Three<*),K3,Dir$,Ct$, Maxval mark , Maxval s< *) 


GOTO  1520 

Four: !  ************************************** 

REDIM  Four<Loc-01dloc,5) 

FOR  1=01 dJ oc  TO  Loc-1 
K4*K4+ 1 
FOR  J=1  TO  5 
Four<K4, J)=Stress<I,J) 

NEXT  J 
NEXT  I 
Loc  =01 dl oc 
FOR  K=1  TO  4 

MAT  SEARCH  Four <*, K+l ), MAX; Maxval s<4, K) 

MAT  SEARCH  Four < * , K  + 1 > , M I N ; Maxval s < 4 , K+4 > 

NEXT  K 

Maxval mark  =4 

CALL  Di  rect  i  onpl  ot<Elem<*>,  Four (*> ,K4,Dir$,Ct$,Maxvalmark, Maxval s( *)  > 
GOTO  1520 

Five: !  ************************************* 

REDIM  Fi ve<Loc-01 dl oc , 5> 

FOR  1=01 dl oc  TO  Loc-1 
K5=K5+ 1 
FOR  J=1  TO  5 
Five<K5,J)=Stress<I,J) 

NEXT  J 
NEXT  I 
Loc=01 dl oc 
FOR  K=1  TO  4 

MAT  SEARCH  Fi ve<*, K+l >, MAX; Maxval s<5, K> 

MAT  SEARCH  Fi ve<*, K+l) , MIN; Maxval s<5, K+4) 

NEXT  K 

Maxval mark=5 

CALL  Directi onpl ot  <  E 1 em<*) , Fi ve<*> , K5, Di rf , Ct  $, Maxval mark , Maxval s ( * ) > 
GOTO  1520 


************************************** 


REDIM  Six(Loc-01d1oc,5> 

FOR  1=01 d 1 oc  TO  Loc-1 
K6=K6+ 1 
FOR  J=1  TO  5 
Six(K6,J)=Stress<I,J) 

NEXT  J 
NEXT  I 
Loc*01d!oc 
FOR  K=1  TO  4 

MAT  SEARCH  Si x< *, K+ 1 ), MAX; Maxval s<6 , K ) 

MAT  SEARCH  S i x< *, K+ 1 ), MIN; Maxval s < 6, K+4 > 

NEXT  K 

Maxval mark  =  6 

CALL  D i rec  t i onp 1 ot  <E 1 em(*> , Si x<») , K6, D i r$, Ct $, Maxval mark , Maxval s<  *  >  > 

GOTO  1520 

Seven: !  *********************************** 

REDIM  Seven<Loc-01 dl oc , 5) 

FOR  1=01 dl oc  TO  Loc-1 

K7=K7+ 1 

FOR  J=1  TO  5 
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2566 

2576 

2586 

2596 

2660 

2610 

2620 

2630 

2640 

2650 

) 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2818 

) 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 


Seven <K7, J)=Stress<I,J> 

NEXT  J 
NEXT  I 
Loc=01 dl oc 
FOR  K»1  TO  4 

MAT  SEARCH  Seven<*,K+l > ,  MAX;  Maxval  s<7,  K) 

MAT  SEARCH  Seven< *, K+ 1 > , M IN; Maxval s < 7, K+4 > 

NEXT  K 

Maxval mark=7 

CALL  Directi  onpl  ot<Elem<*),Seven<*>,K7,  Dirt,  Ct  $,  Max  val mark ,  Maxval  s<*> 
GOTO  1520 

Eight: !  ************************************* 

REDIM  Ei ght <Loc-01 dl oc , 5) 

FOR  1*01 dl oc  TO  Loc-1 

K8=K8+ 1 

FOR  J=1  TO  5 

Eight ( K8 , J ) =St  ress  < I ,  J) 

NEXT  J 
NEXT  I 
Loc  *01 dl oc 
FOR  K= 1  TO  4 

MAT  SEARCH  E i ght < *, K+l > , MAX; Maxval s < 8 , K > 

MAT  SEARCH  E i ght < * , K+ 1 ) , M I N ; Maxval s < 8 , K+4 > 

NEXT  K 

Maxval mark=8 

CALL  Di rec t i onp 1  ot  <E1  em<*> ,  Ei ght  <* > , K8, Di r$, Ct  $, Maxval mark , Maxval s<*) 
GOTO  1520 

Nine: !  ************************************** 

REDIM  Ni ne<Loc-01 dl oc , 5) 

FOR  I ®0 1 d 1 oc  TO  Loc-1 

K9=K9+ 1 

FOR  J=1  TO  5 

Ni ne<K9, J)=Stress< I , J) 

NEXT  J 
NEXT  I 
Loc*01 dl oc 
FOR  K= 1  TO  4 

MAT  SEARCH  N i ne < * , K  + 1 ) , MAX ; Maxval s < 9 , K > 

MAT  SEARCH  Ni ne<*, K+l ) , MIN; Maxval s<9, K+4) 

NEXT  K 

Maxval mark =9 

CALL  Di rec t i onpl ot<Elem<*),Nine<*),K9, Dirt, Ctt, Max valroark, Max vals<*>) 

GOTO  1520 

Ten: !  ************************************** 

REDIM  Ten<Loc-01dloc,5> 

FOR  1*01 dl oc  TO  Loc-1 
K10*K10+1 
FOR  J*1  TO  5 
Ten<K10, J)=Stress( I , J) 

NEXT  J 
NEXT  I 

FOR  K*1  TO  4 

MAT  SEARCH  Ten<*, K+l > , MAX; Maxval s<10, K) 

MAT  SEARCH  Ten< *, K+l >, MIN; Maxval s < 10, K+4) 

NEXT  K 

Maxval mark  *  10 

CALL  Di rec t i onp 1ot<Elem<*>, Ten <*>,K10, Dirt, Ctt, Maxval mark , Maxval s < * ) ) 

GOTO  1520 

CALL  Chart 2< Maxval s<*), Dirt, Ctt, Char tflag) 

MAT  SORT  Stress<*,l> 

!  ************************************************************************* 
INPUT  "More  Data?<Y/N> " .Huht 
IF  Huh$«"N"  THEN  GOTO  3270 
PRINT  PAGE 


40 


3200  PRINT  “Stress  Component  or  Output  Fil  e?(SC/QP) 11 

3210  INPUT  fin* 

3220  IF  <  An*<  >"SC" >  AND  <An*O"0P">  THEN  PRINT  "INVALID  ENTRY" 

3230  IF  <An*<>"SC">  AND  <fin*< >"OP" >  THEN  GOTO  3200 

3240  PRINT  PAGE 

3250  IF  An*® "SC"  THEN  GOTO  Scalemenu 

3260  IF  An*-"OP"  THEN  GOTO  330 

3270  SUBEXIT 

3280  SUBEND 

3290  SUB  Chart 

3300  PLOTTER  IS  "GRAPHICS" 

3310  GRAPHICS 

3320  FOR  1-18  TO  94  STEP  4 

3330  MOVE  1,99 

3340  AREA  COLOR  6/6 ,  (I - 1 8 > -'76 ,  1 

3350  RECTANGLE  4,1, FILL 

3360  NEXT  I 

3370  FOR  1*18  TO  94  STEP  4 

3380  MOVE  1,98 

3390  AREA  COLOR  4/6 , < I - 1 8 ) /76 , 1 

3400  RECTANGLE  4,1, FILL 

3410  NEXT  I 

3420  FOR  1*18  TO  94  STEP  4 
3430  MOVE  1,97 

3440  AREA  COLOR  2/6,  <  I-18W6,  1 

3450  RECTANGLE  4,. 8, FILL 

3460  NEXT  I 

3470  MOVE  0,98 

3480  PEN  1 

3490  CSIZE  2.5, .6 

3500  LABEL  "ZERO  STRESS" 

3510  MOVE  100,98 

3520  LABEL  "YIELD  STRESS" 

3530  CSIZE  2.5,1 

3540  LORG  5 

3550  MOVE  61,94 

3560  LABEL  "EFFECTIVE  STRESS" 

3570  LORG  l 

3580  PRINT  PAGE 

3590  SUBEXIT 

3600  SUBEND 

3610  SUB  Ful 1 <Zmax,A,B,D> 

3620  A*Zmin 

3630  B*Zmax 

3640  D*Zmax/1.23 

3650  SUBEXIT 

3660  SUBEND 

3670  SUB  One_hal f <Zmax, A, B, D> 

3680  A*Zmin 

3690  B-Zmax/2 

3700  D*Zmax/1.23/2 

3710  SUBEXIT 

3720  SUBEND 

3730  SUB  Two_hal f <Zmax, A, B, D) 

3740  A*Zmax/2 

3750  B*Zmax 

3760  D*Zmax/l . 23/2 

3770  SUBEXIT 

3780  SUBEND 

3790  SUB  One_f our t h < Zmax , A, B , D ) 

3800  A*Zmin 

3810  B*Zmax/4 

3820  D-Zmax/1 . 23/4 

3830  SUBEXIT 

3840  SUBEND 

3850  SUB  Two_fourth(Zmax, A, B, D) 
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3660 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 


A*Zmax/4 
B*Zm*x/2 
D«Zmax/l . 23/4 
SUBEXIT 
SUBEND 

SUB  Three_fourth<Zmax, A, B, D) 

A«Zmax/2 
B*3*Zmax/4 
D«Zmax/l . 23/4 
SUBEXIT 
SUBEND 

SUB  Four_fourth<Zmax, A, B, D) 

A«3*Zmax/4 

B=2max 

D*Zmax/l . 23/4 

SUBEXIT 

SUBEND 

SUB  Rscal e<Rmax, A, B, D) 

A*0 

B=Rmax*l . 23*1 . 2 
D*Rmax*l . 2 
SUBEXIT 
SUBEND 

SUB  Di  rec  t  ionplot<Elem<*),Array<*),Kount,Dir$,Ct$,Maxvalmark,Maxuals<*)> 
OPTION  BASE  1 
DIM  Element <50, 3) 

IF  Di rf ="R"  THEN  Dir=l 
IF  Di rt*"Z"  THEN  Dir=2 
IF  Di r**"T"  THEN  Dir=3 
IF  Di r$*"R2"  THEN  Dir*4 


IF  Di rt*"Z"  THEN  Dir=2 
IF  Di r**"T"  THEN  Dir=3 
IF  Di r$*"R2"  THEN  Dir=4 
01 ddi r=Di r 

IF  Ct**"-"  THEN  Dir=Dir+4 
FOR  1*1  TO  Kount 
Elnum*Array< 1,1) 

Hue*El em<E1 num, 11) 

IF  Hue*l  THEN  4430 

IF  Hue>6  THEN  Hue=Hue-4.5 

Hue»Hue/6 

IF  Maxvals<Maxualmark,Dir)*0  THEN  GOTO  4290 

Int ens i t y=Array < 1 , 01 dd i r  +  1 )/Maxval s  <  Max ual mark , Di r  > 

IF  <Array<I,Olddir+l)<0>  AND  <Ct$*"+")  THEN  Intensity*© 

IF  < Array < I,01ddir*l)>0)  AND  CCt**"-")  THEN  Intensity*© 

GOTO  4300 

Intensi ty*0 

PLOTTER  IS  E 1 ement  < * ) 

Luminosi ty*l 

IF  Intensity*©  THEN  Lumi nos i t y= . 15 
IF  Intensity*©  THEN  Intensi ty*. 90 
MOVE  El em<El num, 3) , El em<El num, 4) 

PEN  1 

DRAW  El em<E1 num, 5) , El em<El num, 6) 

DRAW  El em<E 1 num, 7),Elem(EJnum,8) 

DRAW  El em<E 1 num, 9 > , E 1 em<E 1 num, 10) 

DRAW  E 1 em<E 1 num, 3) , E 1 em(E 1 num, 4 ) 

PLOTTER  ElementC*)  IS  OFF 

AREA  COLOR  Hue, Int ens i t y, Lumi nos i t y 

MAT  PLOT  Element, FILL 

NEXT  I 

SUBEXIT 

SUBEND 

SUB  Chart2<Maxuals(*),Dir$,Ctf,Chartflag) 

OPTION  BASE  1 

PRINT  PAGE 

SCALE  0, 123,0, 100 

IF  Di r#*"R"  THEN  Index=l 

IF  Dir#«"Z"  THEN  Index*2 
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IF  Dir*®“T"  THEN  Index=3 

IF  Di r$*"RZ"  THEN  Index=4 

IF  Ct*="-"  THEN  Index® Index+4 

MOVE  €2,94 

CSIZE  2.5, 1 

LORG  5 

IF  <Ct*«"-">  AND  (Di r$="R" )  THEN  LABEL  "COMPRESSIVE  RADIAL  STRESS 

IF  <Ct*®"-">  AND  <Dir*="Z">  THEN  LABEL  "COMPRESSIVE  AXIAL  STRESS" 

IF  <Ct*«"-">  AND  <Dir*®“T“>  THEN  LABEL  "COMPRESSIVE  HOOP  STRESS" 

IF  <Ct*»"-">  AND  <Dir$="RZ">  THEN  LABEL  "NEGATIVE  SHEAR  STRESS" 

IF  <Ct*®"+")  AND  <Dir*="R">  THEN  LABEL  "TENSILE  RADIAL  STRESS" 

IF  <Ct*="+")  AND  <Di r*="Z” )  THEN  LABEL  "TENSILE  AXIAL  STRESS" 

IF  <Ct*®"+">  AND  <Dir*="T">  THEN  LABEL  "TENSILE  HOOP  STRESS" 

IF  <Ct«»"+")  AND  < D i rf  = "RZ" >  THEN  LABEL  "POSITIVE  SHEAR  STRESS" 
LORG  1 

FOR  I®1  TO  10 

Maxval =Maxual s  < I , Index) 

IF  <Ct*®"-")  AND  < Maxval >0 )  THEN  Maxval=0 
IF  <Ct*="+">  AND  (Maxval <0)  THEN  Maxua1=0 
Hue*I 

IF  Hue>6  THEN  Hue®Hue-4.5 
Hue®Huex6 

MOVE  <I-1>*12.4,98 
AREA  COLOR  Hue, 1, 1 
RECTANGLE  2, 2, FILL 
MOVE  <I-1>*12.4*2.2,98 
PEN  1 

CSIZE  2.8, 1x2 
LABEL  Maxual 
NEXT  I 
Chart f 1 ag=l 
SUBEXIT 
SUBEND 


1  !  ********  PROGRAM  "PREPS”  OF  THE  PREPPY  PACKAGE  ******** 

2  !  THIS  IS  A  BLACK  &  WHITE  VERSION  OF  PREP1 

10  OPTION  BASE  1 

20  PRINTER  IS  16 

30  PLOTTER  IS  "GRAPHICS" 

40  GRAPHICS 

50  CALL  Plots 

60  END 

70  SUB  Plots 

80  OPTION  BASE  1 

90  INPUT  “Enter  name  of  Superqrid  f i 1 e: " , Sf i 1 e* 

100  DISP  "*  Reading  in  Supergrid  File:  ";Sfile* 

110  Sfi  le*=Sfi  1e*8<":C12“ 

120  ASSIGN  ttl  TO  Sf i le* 

130  DIM  Elem<800, 11) 

140  READ  #l;Elemax 

150  REDIM  ElemCEIemax, 11) 

160  MAT  READ  #i;Elem 

170  MAT  SEARCH  El em(*, 7) , MAX; Zmaxl 

180  MAT  SEARCH  E 1 em< * , 9 ) , MAX ; Zmax2 

190  MAT  SEARCH  El em<*, 8) , MAX; Rmaxl 

200  MAT  SEARCH  El em<*, 10 > , MAX; Rmax2 
210  Zmax=MAX(Zmaxl , Zmax2> 

220  Rmax=MAX(Rmaxl , Rmax2) 

230  DIM  Stress*<800)C80] 

240  DIM  E1ement(10,3>,Stress<800,7) 

250  REDIM  StressCEl emax, 7> 

260  INPUT  "Material  plot  or  Stress  pi ot?(MxS> " , Choi ce* 

270  IF  Cho i  ce*= "M"  THEN  GOTO  Scalemenu 

280  INPUT  "Enter  name  of  output  file:", Out* 

290  0ut*=0ut*8<" :  C12" 

300  ASSIGN  #3  TO  Out* 

310  DISP  "  *  Organizing  Stress  Array  *" 

320  READ  #3; Stress*< 1 > 

330  FOR  1=1  TO  Elemax 

340  READ  #3; Stress*< I ) 

350  Stressd,  1  )=VAL(St  ress*(I  >  C  1 , 4  3  > 

360  Stress<I,2>=VAL(Stress*<IH5,  1 1 3  > 

370  StressCI, 3)=VAL<Stress*< I >C 1 3, 19 J) 

380  Stress<I,4>=VAL<Stress*<I>l2l,27]> 

390  Stress<I,5)=VAL<Stress*<I)C29,35]) 

400  Stress( I , 6) =VAL (Stress* < 1)137,433) 

410  Stress( I , 7>=E1 em( I , 1 1 > 

420  NEXT  I 


430 

Scalemenu:  ! 

440 

PRINT  "Scale  Menu: 

FULL 

H 

1  450 

PRINT  " 

1x2 

II 

460 

PRINT  " 

2x2 

II 

470 

PRINT  " 

1x4 

II 

480 

PRINT  " 

2x4 

II 

490 

PRINT  " 

3x4 

II 

500 

PRINT  " 

4x4 

II 

1  510 

INPUT  Menu* 

1  520 

PRINT  PAGE 

530 

GCLEAR 

548 

IF  <Zmax>Rmax) 

AND 

CMenu*="FULL" 

) 

THEN  CALL 

Ful  1  (Zma  >,  B,  D) 

550 

IF  (Zmax>Rmax) 

AND 

(Menu*-" 1 x2 

•) 

THEN  CALL 

One_halff  <ax,A,B,D> 

560 

IF  <Zmax>Rmax> 

AND 

(Menu*="2x2 

’ ) 

THEN  CALL 

Two_half  .max,A,B,D) 

!  570 

IF  (Zmax>Rmax) 

AND 

(Menu*-" 1x4 

' ) 

THEN  CALL 

One_foi  h(Zmax, A, B, D> 

1  580 

IF  (Zmax>Rmax) 

AND 

<Menu*-"2x4 

' ) 

THEN  CALL 

Two_fo  th(Zmax, A, B, D) 

590 

IF  (Zmax>Rmax) 

AND 

( Menu*- " 3x4 " ) 

THEN  CALL 

Three  JurthCZmax, A, B, D> 

600 

IF  (Zmax>Rmax) 

AND 

(Menu*-"4x4 

' ) 

THEN  CALL 

Four  Vourth(Zmax, A, B, D) 

610 

IF  . 99*Zmax< =Rmax 

THEN  CALL 

Rsc  al e(Rmax,  A,  B 

i  D) 

620 

SCALE 

A,  B, 

0,  D 

630 

IF  Choice* 

« " M "  THEN 

GOTO 

Ef f ec  t i we 

640  PRINT  "Specify  Stress  Component :  R, Z, T, RZ, EFF" 

52 


650  INPUT  Dir* 

660  IF  Chartf 1 ag=l  THEN  SCALE  A,B,0,D 

670  IF  <Di r$<  >"R" )  AND  <Di r*< >"Z" )  AND  < Di r*< > "T " >  AND  <Dir$<>"RZ")  AN 

D  <Dir*<)"EFF")  THEN  PRINT  "INVALID  ENTRY" 

680  IF  <Dir*<>"R")  AND  <Di r*< >"Z" )  AND  <Di r*< >"T"  )  AND  <Dir*<>"R2")  AN 

D  <Dir*<>"EFF")  THEN  OUTPUT  7, 6; "BP7, 300, 5; BP6, 300, 5" 

690  IF  <Dir$<>"R")  AND  <Dir*<>"2">  AND  <Dir*<>"T")  AND  <Di r$< >"RZ" )  AN 

D  <Dir»<>"EFF")  THEN  GOTO  640 

700  IF  D i r*<  > "EFF"  THEN  GOTO  740 

710  PRINT  PAGE 

720  PRINT  "Component  =  EFF" 

730  GOTO  Effective 

740  PRINT  "Specify  Direction: <-/+)" 

750  INPUT  Ct* 

760  PRINT  PAGE 

770  GCLEAR 

780  PRINT  "Component  =  " ; Dir*, "Direct i on  =  ";Ct* 

790  GOTO  Other  direction 


800  Effective 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110  ! 

1120  ! 

1130  ! 

1140 
1150 
1 160 
1170 
1180 
1190 
1200  ! 

1210  ! 

1220  ! 

1230 

1231 

1232 

1233 

1234 


IF  Choi ce*="S"  THEN  CALL  Chart 

SCALE  A, B, 0, D 

MAT  SORT  Stress<*,l) 

FOR  1*1  TO  Elemax 
Hue*El em< 1,11) 

IF  Hue=l  THEN  1190!***IF  YOU  WANT  AIR  PLOTTED  COMMENT  THIS  LINE 

IF  Choi ce*="M"  THEN  GOTO  910 

Intensity=Stress<I,6)/£lem(I,2) 

Intense=Intensity 
IF  Intensity>l  THEN  Intensity=l 
IF  lntensity<0  THEN  Intensity*© 

IF  Hue>6  THEN  Hue=Hue-4.5 
Hue=Hue/6 

PLOTTER  IS  E 1 ement  < * ) 

MOVE  Elem<I,3),E1em<I,4) 

PEN  -1 

DRAW  Elem<I,5),Elem<I,6) 

DRAW  Elem<I,7),Elem<I,8) 

DRAW  Elem<I,9),E1em(I, 10) 

DRAW  Elem(I,3),Elem<I,4) 

PLOTTER  Elements)  IS  OFF 

IF  Choi ce*="M"  THEN  I nt ens i t y =Hue 

AREA  COLOR  1 , 1 , I nt ens i t y 

MAT  PLOT  Element, FILL 

IF  Choi ce*="S"  THEN  1120 

MOVE  Elem(I,3),Elem<I,4) 

PEN  -1 

DRAW  Elem<I,5),Elem<I,6) 

DRAW  Elem(I,7),E!em(I,8) 

DRAW  Elem<I,9),Elem<I, 10) 

DRAW  Elem(I,3),E1em(I,4) 

*****  ELEMENTS  ABOVE  YIELD  ***** 

IF  ( Intense< 1.02)  OR  <Di r*<  >"EFF" )  THEN  GOTO  1190 
MOVE  Elem<I,3),E!em<I,4) 

DRAW  Elem<I,7),Elem<I,8) 

MOVE  Elem<I,5),Elem(I,6) 

DRAW  Elem<I,9),Elem<I, 10) 

NEXT  I 

******************************** 

PRINTER  IS  0 

INPUT  "DUMP  GRAPHICS?  <Y/N)",Dg$ 

IF  Dg*="Y"  THEN  DUMP  GRAPHICS 
IF  Dg*="Y"  THEN  PRINT 
PRINTER  IS  16 


S3 


1240  IF  Choi ce*="M"  THEN  INPUT  “Look  at  another  sect i on?< Y/N> " , Lo 

ok* 


1250  IF  Look*="N"  THEN  3210 

1260  IF  Choice*="M"  THEN  GOTO  Scalemenu 

1270  INPUT  "More  data  <Y/NV,Huh* 

1280  IF  Huh*="N"  THEN  GOTO  3210 

1290  PRINT  PfiGE 

1300  PRINT  "Stress  Component  or  Output  File:  SCxQP" 

1310  INPUT  fin* 

1320  IF  <fin*<>"SC">  FIND  <An*<  >"0P" )  THEN  PRINT  "  INVfiLID  ENTRY" 

1330  IF  <fln*<  >"SC" )  AND  <An*O"0PM>  THEN  GOTO  1300 

1340  PRINT  PAGE 

1350  IF  fln*="SC"  THEN  GOTO  Scalemenu 

1360  IF  An**"0P"  THEN  GOTO  280 

1370  Other__di  rec  t  i  on:  ! 

1380  DIM  One<350,5>,Two<850,5>,Three<300,5>,Four<850,5>,Fi ve<300,5>,Six<850 
5> , Seven <850, 5) , Ei ght  <300, 5),Nine<300,5), Ten <850, 5>,Maxvals<10,8> 

1390  K1=K2=K3=K4=K5=K6=K7=K8=K9=K10=0 

1400  MAT  SORT  Stress<*,7>  DES 

1410  Loc*£1emax+l 

1420  FOR  Mat  1*1  TO  10 

1430  MAT  SEARCH  St ress < * , 7 > , LOC < =Mat 1 ) ; 0 1 d 1 oc 

1440  IF  01 dl oc=El emax+ 1  THEN  1460 

1450  ON  Stress<Loc-l , 7)  GOTO  One,  Two, Three , Four , Fi ve, Si x, Seven, Ei ght , Ni ne, T 

n 


1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

I860 


I  CONTINUE 
NEXT  Mat  1 
GOTO  3080 

One: !  ************************************** 

REDIM  0ne<Loc-01dloc,5) 

FOR  1*01 dl oc  TO  Loc-1 
K 1 *K 1 + 1 
FOR  J*1  TO  5 
One<Kl ,J)=Stress<I,J> 

NEXT  J 
NEXT  I 
Loc=01dloc 
FOR  K* 1  TO  4 

MAT  SEARCH  One < * , K  + 1 > , MAX ; Maxval s < 1 , K ) 

MAT  SEARCH  One <*, K+ 1 >, M I N ; Maxval s < 1 , K+4 > 

NEXT  K 

Maxval mark*l 

CALL  Di  rec  t  i  onpl  ot<Elem<*),One<*>,Kl,Dir*.Ct*, Maxval mark,Maxvals<*)> 
GOTO  1460 

Two:  !  ************************************ 

REDIM  Two<Loc-01dloc,5> 

FOR  1*01 dl oc  TO  Loc-1 

K2*K2+1 

FOR  J*1  TO  5 

Two<K2, J)=Stress< I , J> 

NEXT  J 
NEXT  I 
Loc*01 dl oc 
FOR  K* 1  TO  4 

MAT  SEARCH  Two < * , K  + 1 ) , MAX ; Max val s < 2 , K > 

MAT  SEARCH  Two< *, K+l ), MIN; Maxval s <2, K+4> 

NEXT  K 

Maxval mark  *2 

CALL  Di rec t i onp lot<Elem<*), Two <*>,K2, Dir*, Ct$, Maxval mark ,Maxvals<*>> 

GOTO  1460 

Three: !  f*********************************** 

REDIM  Three<Loc-01dloc,5> 

FOR  1*01 dl oc  TO  Loc-1 
K3-K3+1 
FOR  J*1  TO  5 
Three <K3, J)*Stress< I , J) 


54 


1870 

1888 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

) 

I960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 


NEXT  J 
NEXT  I 
Loc*01 dl  oc 
FOR  K* 1  TO  4 

MAT  SEARCH  Three < *, K  + 1 ) , MAX; Maxval s <3 , K> 

MAT  SEARCH  Three <* , K+ 1 > , MIN; Maxval s  <  3, K+4  > 

NEXT  K 

Maxval mark  =3 

CALL  Di  rec t  i  onpl  ot  <E1  em(*>  ,  Three < * >  ,  K3,  Di  r$ ,  Ct  #,  Maxval  mark  ,  Maxval  s< *  > 
GOTO  1460 

Four: !  ***»*•*•*••**•****•**•*•**#***  **************** 

REDIM  Four(Loc-01dloc,5> 

FOR  I *01 d) oc  TO  Loc-1 

K4*K4+ 1 

FOR  J*1  TO  5 

Four <K4, J)=Stress<I,J) 

NEXT  J 
NEXT  I 
Loc*01 dl oc 
FOR  K* 1  TO  4 

MAT  SEARCH  Four  (*,  K+l  ),  MAX;  Maxval  £<'4,  K> 

MAT  SEARCH  Four < * , K+ 1) , M I N ; Maxval s < 4 , K+4 > 

NEXT  K 

Maxval mark  *4 

CALL  Di  rec  t  i  onpl  ot  <E  1  em( *> ,  Four  < *) ,  K4,  Di  r$,  Ct  $,  Maxval  mark  ,  Maxual  s < * > ) 

GOTO  1460 

Fioe:  !  ************************************* 

REDIM  Fi ueCLoc-01 dl oc , 5) 

FOR  I*01dloc  TO  Loc-1 

K5*K5+1 

FOR  J=1  TO  5 

Five<K5,J)=Stress(I,J) 

NEXT  J 
NEXT  I 
Loc*01 dl oc 
FOR  K-l  TO  4 

MAT  SEARCH  F i ve<*, K+ 1 ), MAX; Maxval s < 5, K > 

MAT  SEARCH  F i ve < * , K+ 1 > , M I N ; Maxual s C 5 , K+4 ) 

NEXT  K 

Maxval mark  =5 

CALL  Di rec  t i onp 1 ot  <E 1 em(*) , Fi ve<*) , K5, Di r$ , Ct  $, Maxval mark , Maxval s < * ) ) 

GOTO  1460 

Six: !  ************************************** 

REDIM  Six<Loc-01d1oc,5> 

FOR  I *01 d 1 oc  TO  Loc-1 
K6*K6+ 1 
FOR  J*1  TO  5 
Six<K6,J>*Stress<I,J> 

NEXT  J 
NEXT  I 
Loc*01 dl oc 
FOR  K* 1  TO  4 

MAT  SEARCH  Si x<*, K+l >, MAX; Maxval s<6, K > 

MAT  SEARCH  Si x<*, K+l ), MIN; Maxval s<6, K+4 > 

NEXT  K 

Maxval mark  *6 

CALL  D i rec t i onpl ot (El em<*> , Si x<*> , K6, Di r$, Ct$, Maxval mark , Maxval s  <  *  >  > 

GOTO  1460 

Seven: !  *********************************** 

REDIM  Seven<Loc-01 dl oc , 5) 

FOR  1*01 dl oc  TO  Loc-1 

K7-K7+1 

FOR  J*1  TO  5 

Seven <K7, J)*Stres*< I , J) 

NEXT  J 


I 


2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

> 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

> 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 


NEXT  I 
Loc=01 dl oc 
FOR  K* 1  TO  4 

MAT  SEARCH  Seven<*, K+l > , MAX; Maxval s<7, K> 

MAT  SEARCH  Seven< * , K+ 1 ) , M I N; Maxval s < 7, K+4 ) 

NEXT  K 

Maxval mark=7 

CALL  Directionp1ot<E1em(*)>Sewen<*),K7,Dir$,Ct$,Maxwalmark>Maxvals(*) 

GOTO  1460 

Eight: !  ************************************* 

REDIM  Ei ght <Loc-01 dl oc , 5) 

FOR  1*01 dl oc  TO  Loc-1 

K8*K8+1 

FOR  J*1  TO  5 

Eight  <K8,J)=Stress<I,J> 

NEXT  J 
NEXT  I 
Loc*01 dl oc 
FOR  K*1  TO  4 

MAT  SEARCH  E i ght < * , K  + 1 > , MAX; Maxval s < 8 , K > 

MAT  SEARCH  E i ght <*, K+ 1 ) , MIN; Maxval s < 8 t K+4 ) 

NEXT  K 

Maxval mark =8 

CALL  Di rect i  onpl  ot  <E1  em<*>  ,  Ei  ght  (*>  ,  K8,  Di  r$ ,  Ct  $,  Maxval  mark  ,  Max val  s  (.* > 
GOTO  1460 

Nine: !  ************************************** 

REDIM  Nine<Loc-01dloc,5> 

FOR  1-01 dl oc  TO  Loc-1 
K9=K9+ 1 
FOR  J*1  TO  5 
Nine<K9, J)=Stress<I, J) 

NEXT  J 
NEXT  I 
Loc=01 dl oc 
FOR  K* 1  TO  4 

MAT  SEARCH  N i ne < * , K+ 1 > , MAX ; Max val s < 9 , K > 

MAT  SEARCH  N i ne < * , K+ 1 ) , M I N; Maxval s < 9 , K+4 > 

NEXT  K 

Maxval mark  =9 

CALL  Di rect i onpl ot<Elem<*),Nine<*),K9,Dir$,Ct*,Maxvalmark, Maxval s<*>  > 
GOTO  1460 

Ten: !  ************************************** 

REDIM  Ten<Loc-01dloc,5) 

FOR  1*01 dl oc  TO  Loc-1 
K10-K10+1 
FOR  J*1  TO  5 
Ten<K10, J)=Stress< I , J) 

NEXT  J 
NEXT  I 

FOR  K* 1  TO  4 

MAT  SEARCH  Ten< *, K+ 1 >, MAX; Maxval s < 1 0, K ) 

MAT  SEARCH  Ten< *, K+l ), MIN; Maxval s < 10 , K+4 > 

NEXT  K 

Maxval mark  *  1 8 

CALL  Directionp1ot<Elefn<*),Ten(*>,K10,Dir$,Ct$,Maxvalmark,Maxvals(*>) 

GOTO  1460 

CALL  Chart  2 < Maxval s<*>,Dirf,Ct$, Chart  flag) 

MAT  SORT  Stress<«, 1 > 

!  ************************************************************************* 
INPUT  "More  Data?<Y/N) " . Hvh« 

IF  Huh$*"N"  THEN  GOTO  3210 
PRINT  PAGE 

PRINT  "Stress  Component  or  Output  F i 1 e?< SC/QP ) " 

INPUT  An$ 


56 


K 


fi 


R 


3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 

3650 

3660 

3670 

3680 

3690 

3700 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 


IF  (An*<  >"SC"  >  AND  (An*<  VOP"  >  THEN  PRINT  "INVALID  ENTRY" 
IF  (An*<  >"SC" >  AND  (An*O"0P">  THEN  GOTO  3140 

PRINT  PAGE 

IF  An*="SC"  THEN  GOTO  Scalemenu 
IF  An** "OP"  THEN  GOTO  280 
SUBEXIT 
SUBEND 
SUB  Chart 

PLOTTER  IS  "GRAPHICS" 

GRAPHICS 

FOR  1*18  TO  94  STEP  4 
MOVE  1,97 

AREA  COLOR  1,1, (1-18) '76 
RECTANGLE  4, 4, FILL 
NEXT  I 
MOVE  0,98 
PEN  1 

CSIZE  2.5, .6 
LABEL  "ZERO  STRESS" 

MOVE  100,98 

LABEL  “YIELD  STRESS" 

CSIZE  2.5 
LORG  5 
MOVE  61,94 

LABEL  "EFFECTIVE  STRESS  X  >  2>.  ABOVE  YIELD" 

LORG  1 
PRINT  PAGE 
SUBEXIT 
SUBEND 

SUB  Ful 1 (Zmax, A, B, D) 

A-0 

B*Zmax 
D=Zmax/ 1 . 23 
SUBEXIT 
SUBEND 

SUB  One_hal f (Zmax, A, B, D> 

A=0 

B*Zmax/2 
D*Zmax/ 1 . 23x2 
SUBEXIT 
SUBEND 

SUB  Two_hal f (Zmax, A, B, D) 

A*Zmax/2 

B*Zmax 

D=Zmax/l . 23/2 

SUBEXIT 

SUBEND 

SUB  One_fourth(Zmax, A, B, D) 

A=0 

B*Zmax/4 
D*Zmax/l . 23/4 
SUBEXIT 
SUBEND 

SUB  Two_fourth(Zmax, A, B, D> 

A*Zmax/4 
B*Zmax/2 
D*Zmax/ 1 . 23/4 
SUBEXIT 
SUBEND 

SUB  Three_four th(Zmax, A, B, D> 

A*Zmax/2 
B*3*Zmax/4 
D*Zmax/ 1 . 23/4 
SUBEXIT 
SUBEND 

SUB  Four  fourth(Zmax, A, B, D) 

■S7 


3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 

3 

4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

4420 


ft 


4430 

4440 


ft 


4430 

4460 


ft 


A*3*Zmax/4 

B*Zmax 

D*Zmax/ 1 . 23^4 

SUBEXIT 

SUBEND 

SUB  Rsc al eC Rmax , ft, B, D) 

A*0 

B*Rmax*l . 23*1 . 2 
D*Rmax* 1 . 2 
SUBEXIT 
SUBEND 

SUB  Di  rec t  i  onp  1  ot  CE1  emC*)  ,  Array C*)  ,  (Count ,  Di  r$ ,  Ct$,  Maxval  mark  ,  Maxval  sC * )  > 

OPTION  BASE  1 
DIM  El ement C 10, 3) 

IF  Di r$*"R"  THEN  Dir*l 
IF  Di  r$*"Z“  THEN  Dir  =  2 
IF  Dir*="T"  THEN  Dir  =  3 
IF  Di  r$*"RZ"  THEN  Dir  =  4 
01 ddi r*Di r 

IF  Ct»="-"  THEN  Dir=Dir+4 
FOR  1*1  TO  (Count 
El num*ArrayC I , 1 ) 

Hue=E1 emCEl num, 11) 

IF  Hue* 1  THEN  4250 

IF  Hue>6  THEN  Hue*Hue-4.5 

Hue=Hue/6 

IF  Maxual s< Maxual mark , D i r ) =0  THEN  GOTO  4130 
Intensity=flrray<I,01ddir+l)/Maxoals<Maxoalmark,Dir> 

IF  <flrray<I,Olddir+l)<0)  AND  <Ct*=,,  +  ,,>  THEN  Intensity*© 

IF  <Array(I,Olddir+l)>0>  AND  <Ct*="-")  THEN  Intensity*© 

GOTO  4140 

Intensi ty*0 

PLOTTER  IS  ElementC*) 

MOVE  El emCEl num, 3) , El emCEl num, 4) 

PEN  1 

DRAW  El emCEl num, 5) , El emCEl num, 6) 

DRAW  El emC E 1 num , 7 ) , E 1 em CE 1 num, 8) 

DRAW  ElemCElnum,9),ElemCElnum, 10) 

DRAW  El emCE 1 num, 3 ) , E 1 emC E 1 num, 4 ) 

PLOTTER  ElementC*)  IS  OFF 

IF  I nt ens i t y< . 08  THEN  Intensi ty=. 08 

AREA  COLOR  1 ,  1 , I nt ens i t y 

MAT  PLOT  Element, FILL 

NEXT  I 

SUBEXIT 

SUBEND 

SUB  Chart 2C Maxual sC*), Dirt, Ct$, Char tflag) 

OPTION  BASE  1 

PRINT  PAGE 

SCALE  0, 123,0, 100 

IF  Di r**"R"  THEN  Index*l 

IF  Di r**"Z"  THEN  Index=2 

IF  D i r  $  * " T "  THEN  Index*3 

IF  Dirts'^Z"  THEN  Index=4 

IF  Ctf*"-M  THEN  Index=Index+4 

MOVE  62,98.75 

CSIZE  2.5,1 

LORG  5 

IF  CCt$*"-H)  AND  CDi r$*"R" )  THEN  LABEL  "COMPRESSIVE  RADIAL  STRESS" 

IF  CCt$*"-")  AND  CDi r$*"Z" )  THEN  LABEL  "COMPRESSIVE  AXIAL  STRESS" 

IF  CCt**"-")  AND  CDi r$="T" )  THEN  LABEL  "COMPRESSIVE  HOOP  STRESS" 

IF  CCt**"-")  AND  CDi r$="RZ" )  THEN  LABEL  "NEGATIVE  SHEAR  STRESS" 

IF  CCt**"+")  AND  CDi r$="R" )  THEN  LABEL  "TENSILE  RADIAL  STRESS" 

IF  CCt**"+")  AND  CDi r$*"Z" )  THEN  LABEL  "TENSILE  AXIAL  STRESS" 

IF  CCtf-"*")  AND  CDi r$="T" )  THEN  LABEL  "TENSILE  HOOP  STRESS" 

IF  CCt*-"+")  AND  CDir*="RZ")  THEN  LABEL  "POSITIVE  SHEAR  STRESS" 


MAT  #2  MAT  #3  MAT  #4  MAT  #5  MAT  #6  MAT  #7  MAT  #8 


THEN  Maxval=0 
THEN  Maxval=0 


LORG  1 
MOVE  0,95 
CSI2E  2.525 

LABEL  "MAT  #1  MAT  #2  MAT  #3  M 
MAT  #10" 

FOR  1*0  TO  9 

Maxval =Maxva1 s< 1  +  1 ,  Index) 

IF  <Ct#*"-")  AND  (Maxval >0)  THEN 
IF  <Ct#="+">  AND  <Maxva1 <0>  THEN 
MOVE  1*12,92 
PEN  1 

CSIZE  2.8, 1/2 
LABEL  Maxval 
NEXT  I 

INPUT  “DUMP  GRAPHICS?  <Y/N>",Dg$ 

PRINTER  IS  0 

IF  Dg$*“Y"  THEN  PRINT 

IF  Dg*="Y"  THEN  DUMP  GRAPHICS 

PRINTER  IS  16 

Chartf  1  ag*l 

SUBEXIT 

SUBEND 


!  ********  PROGRAM  "GRID"  OF  THE  PREPPY  PACKAGE  ******** 
OPTION  BASE  1 
PRINTER  IS  16 

INPUT  "Enter  name  of  Grid  file: ".File* 

Name**F i  1  et 
Fi  1e**Fi  1e*8.":C12" 

DIM  P 1 ot  <2000, 3) 

INPUT  "Enter  Devi ce: CRT<C>  or  PI  otter <P> ",  Dev* 

PLOTTER  IS  "GRAPHICS" 

GRAPHICS 

IF  Dev*»"C"  THEN  GOTO  190 

INPUT  "Report<R>  Mini(M)  Oversize  <0>",Size* 

IF  Si ze**"R"  THEN  OUTPUT  7, 5; " IP1000, 0, 9000, 646?" 

IF  Si ze**"M"  THEN  OUTPUT  7, 5; " IP3000, 4000, 5000, 5617" 

IF  Si ze**"0"  THEN  OUTPUT  7, 5; " IP1000, 0, 15102, 1 1400" 

PLOTTER  IS  "9872A" 

IF  Si  ze*="R"  THEN  INPUT  "  Vi  eu>qraph<V)  or  Paper  <P)  ",  Dp* 

OUTPUT  7,5;"VS15" 

IF  Dp**"V"  THEN  OUTPUT  7,5;"VS5“ 

ASSIGN  «1  TO  File* 

MAT  READ  ttl;Plot 

MAT  SEARCH  PI Ot <*, 1) , MAX; Max 

MAT  SEARCH  PI ot <*, 1 > , MIN; Mi n 

MAT  SEARCH  PI ot <*, 2) , MAX; Mar 

S*<  l-Mar/Max>/2 

IF  Mar >Max  THEN  Max=Mar*1.2 

Max*Max-Min 

GCLEAR 

INPUT  "PLOT  MENU:  FULL,  1/2,  2/2,  1/4,  2/4,  3/4,  4/4", Menu* 
IF  Menu**" FULL"  THEN  290 
S*S/VAL<Menu*C33  > 

IF  Menu*="FULL"  THEN  SCALE  Mi n, Max, -S*Max, Cl-S>*Max 

IF  Menu**" 1/2"  THEN  SCALE  Mi n, Max/2, -S*Max/2, < l-S>*Max/2 

IF  Menu*»"2/2"  THEN  SCALE  Max/2,  Max, ~S*Max/2, < l-S>*Max/2 

IF  Menu**" 1/4"  THEN  SCALE  Mi n, Max/4, -S*Max/4, < l-S>*Max/4 

IF  Menu**"2/4"  THEN  SCALE  Max/4 , Max/2 , -S*Max/4 , C 1 -S > *Max/4 

IF  Menu**"3/4"  THEN  SCALE  Max/2, 3*Max/4, -S*Max/4, (l-S>*Max/4 

IF  Menu** "4/4"  THEN  SCALE  3*Max/4, Max, -S*Max/4, < l-S)*Max/4 

MAT  PLOT  Plot 

PRINT  PAGE 

SCALE  0, 123,0, 100 

MOVE  62,2 

CSIZE  3.5,1 

LORG  5 

LABEL  "Grid  File  *  " ;Name* 

W-1.5 

BEEP 

WAIT  300/M 
BEEP 

WAIT  150/M 
BEEP 

WAIT  150/M 
BEEP 

WAIT  300/M 
BEEP 

WAIT  600/M 
BEEP 

WAIT  300/M 
BEEP 
PEN  0 

INPUT  "Look  at  another  sect i on?< Y/N) " , L* 

IF  L»*"Y"  THEN  GOTO  250 
END 


,  »V_  .4  ■'a  »V1  (*t  ■  ««  «•.  t%.  df,  »  «  »*,.  f„  ■  .  |»,  I 


!  ********  PROGRAM  "ELNUM"  OF  THE  PREPPY  PACKAGE  ******** 
OPTION  BASE  1 
PRINTER  IS  16 
PLOTTER  IS  “GRAPHICS" 

GRAPHICS 
CALL  Plots 
END 

SUB  Plots 
OPTION  BASE  1 

INPUT  "Enter  name  of  Superqrid  f i 1 e : “ , Sf i 1 e$ 

Sf i 1 e**Sf i 1 e$&" : C12" 

ASSIGN  *1  TO  Sf i le* 

DIM  Elem(800, 11) 

READ  #1 ; El emax 
REDIM  Elem(Elemax, 11) 

MAT  READ  ttlJElem 

MAT  SEARCH  El em(*, 7) , MAX; Zmax 

MAT  SEARCH  E 1 em ( * , 8 ) , MAX ; Rmax 

PRINT  “Prepare  platen  for  Large-Size  paper,  then  CONT" 
PAUSE 

INPUT  "Plot  Menu:  Full,  1/2,  2/2,  1/4,  2/4,  3/4,  4/4", Mem 


then  CONT" 


3/4,  4/4", Menu* 


(Menu*="FULL" )  THEN  CALL  Ful 1 (Zmax, A, B, D> 
(Menu**11 1 /2" )  THEN  CALL  One_hal f (2max, A,  B,  D) 
<Menu*="2/2">  THEN  CALL  Two  hal f (Zmax, A, B, D) 


<Menu*="2/2"> 
<Menu*="l/4"> 
(Menu*="2/4"> 
<Menu*=“3/4" ) 
(Menu*="4/4"  ) 


IF  (Zmax>Rmax>  AND  (Menu*=“FULL" >  THEN  CALL  Fu 
IF  (Zmax>Rmax>  AND  (Menu**" 1 /2" )  THEN  CALL  One 

IF  (Zmax)Rmax)  AND  (Menu$="2/2" >  THEN  CALL  Two 

IF  (Zmax)Rmax)  AND  (Menu**" 1 /4" >  THEN  CALL  One 

IF  <Zmax>Rmax)  AND  (Menu*="2/4" >  THEN  CALL  Two 

IF  (Zmax>Rmax>  AND  (Menu*="3/4" )  THEN  CALL  Thr 

IF  (Zmax>Rmax>  AND  (Menu*="4/4" >  THEN  CALL  Fou 

IF  . 99*Zmax<=Rmax  THEN  CALL  Rscal e(Rmax, A, B, D> 
SCALE  A, B, 0, D 

OUTPUT  7,5*,  “  I P 1 000 , 0 ,  15102,  11400" 
PLOTTER  IS  " 9872A" 

SCALE  A, B, 0, D 

FOR  1*1  TO  Elemax 

MOVE  El  em( 1 , 3)+E, E 1 em( 1 , 4>+E 

PEN  1 

DRAW  Elem(I,5)+E,Elem(I,6)-E 


THEN  CALL  Two_hal f (Zmax, A, B, D) 

THEN  CALL  One_f ourth(Zmax, A, B, D) 
THEN  CALL  Two_fourth(Zmax, A, B, D) 
THEN  CALL  Three_f ourth(Zmax, A, B, D) 


THEN  CALL  Four  fourthCZmax, A, B, D) 


i<  1 , 4>+E 


DRAW  Elem(I,5)+E,Elem(I,6)-E 
DEG 

DRAW  Elem(I,7)-E,Elem(I,8)-E 
DRAW  Elem(I,9)-E,Elem(I, 10>+E 
DRAW  Elem(I,3)+E,Elem(I,4)+E 
Ldi r=0 
LORG  5 

IF  E1em(I,6)-Elem(I,4)>Elem(I,7)-Elem(I,5>  THEN  Ldir=90 
LDIR  Ldi r 
CSIZE  1 

Jim*<Elem<I,6)-E1em<I,4))/<B-A)*123 
Jer*<El em< I , 7)-El em< I , 5) >/D*100 

IF  (Ldi r*0)  AND  ((Jim<.85)  OR  (Jer<2.3>>  THEN  540 

IF  (Ldi r=90>  AND  ((Jer<.85>  OR  (Jim<2.3>)  THEN  540 

Yax*(MAX(El em(I,6),Elem(I,8))+MIN(E1em(I,4),Elem(I, 10)))/2 

Xax*(MAX(El  em( 1 , 7) , El em( 1 , 9) )+MIN(El em( 1 , 3) , El em( 1 , 5) ) >/2 

MOVE  Xax, Yax 

LABEL  I 

NEXT  I 

SUBEXIT 

SUBEND 

SUB  Ful 1 (Zmax, A, B, D) 

A-0 

B*Zmax 
D*Zmax/l . 23 
SUBEXIT 
SUBEND 

SUB  One  hal f (Zmax, A, B, D) 

A*0 

B«Zmax/2 
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D«Zmax/l . 23/2 

SUBEXIT 

SUBEND 

SUB  Two_hal f <Zmax, fl, B, D> 

ft«Zaax/2 

B»Zaax 

D*Zmax/l . 23/2 

SUBEXIT 

SUBEND 

SUB  On*  fourth (Zmax, R, B, D) 

fl«0 

BaZmax/4 
D>Zm*x/l . 23/4 
SUBEXIT 
SUBEND 

SUB  Two_fourth<Zmax, fl, B, D) 

fi**Zmax/4 

B*Zmax/2 

D«Zmax/l . 23/4 

SUBEXIT 

SUBEND 

SUB  Thr**_fourth<Zmax,fl,  B,D> 

R»Zmax/2 

B*3*Znax/4 

D>Zmax/ 1.23/4 

SUBEXIT 

SUBEND 

SUB  Four_f ourth<Zmax,  A,  B,  D> 

9«3*Z*ax/4 

BaZmax 

D>Zmax/l . 23/4 

SUBEXIT 

SUBEND 

SUB  Rscal *<Rmax, fl, B,  D) 

n>0 

B-Rmax* 1.23* 1.2 
D*Rmax*l . 2 
SUBEXIT 
SUBEND 
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